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The stress-energy tensor of a quark moving through a strongly coupled A/" = 4 supersymmetric 
Yang-Mills plasma is evaluated using gauge/string duality. The accuracy with which the resulting 
wake, in position space, is reproduced by hydrodynamics is examined. Remarkable agreement is 
found between hydrodynamics and the complete result down to distances less than 2/T away from 
the quark. In performing the gravitational analysis, we use a relatively simple formulation of the 
bulk to boundary problem in which the linearized Einstein field equations are fully decoupled. Our 
analysis easily generalizes to other sources in the bulk. 



I. INTRODUCTION 

The discovery that the quark-gluon plasma produced 
in heavy ion collisions at RHIC behaves as a nearly ideal 
fluid [TJ E] has prompted much interest in understanding 
the dynamics of strongly coupled non-Abelian plasmas. 
Much recent theoretical work has explored the dynamics 
of maximally supersymmetric Yang-Mills (A/" = 4 SYM) 
plasma. (See, for example, Refs. [SimiHEllIllHlEllini 
[TUIHIIIIITIITSIIISII^.) Thanks to gauge/string (or 
AdS/CFT) duahty [H [19], the properties of a strongly 
coupled A/" = 4 SU{Nc) SYM plasma, in the limit of large 
iVc, are under much better theoretical control than is the 
case for a strongly coupled QCD plasma. Available evi- 
dence from the equation of state [TS] , screening lengths 
[12], and viscosity |6l |2Q], suggests that a strongly cou- 
pled, large- iVc, Af = 4 SYM plasma mimics the properties 
of a real QCD plasma in the temperature range relevant 
to heavy ion collisions (roughly 1.5 to 4 times the de- 
confinement temperature) sufficiently well that one may 
view JV — 4: SYM plasma as a useful model system for 
QCD. 

Heavy quarks produced by hard processes during the 
early stages of a heavy ion collision may transfer much 
of their energy and momentum to the medium if they 
travel a sufficient distance through the fireball before po- 
tentially escaping (and hadronizing) . Analysis of the dis- 
tribution and correlations in produced jets can provide 
information about the dynamics of the plasma, includ- 
ing the rates of energy loss and momentum broadening 
of quarks traversing the plasma [^TJ [2U [53] • 

A heavy quark moving through a plasma will disturb 
the surrounding medium and its motion will, in turn, be 
influenced by the medium. As the quark moves, frictional 
forces will transfer energy and momentum from the quark 
to the plasma. A natural question to consider is where 
does the energy and momentum lost by the quark go? In 
other words, what is the form of the wake, as deflned by 
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the change in the expectation value of the stress-energy 
tensor, 

AT^-(a;) = (r'^''(x))with quark- (T^'^(a:))w/o quark- (1.1) 

For distances asymptotically far from the quark, one may 
address this question using hydrodynamic approxima- 
tions. This approach was used in Ref. [H]. However, 
in this work the form of the effective sources for hydro- 
dynamics were not fully specified. 

Gauge/string duality allows one to compute many ob- 
servables probing non-equilibrium dynamics of strongly 
coupled A/" = 4 supersymmetric Yang-Mills theory, in- 
cluding the rate of energy loss of a heavy quark moving 
through an SYM plasma [4]. (See also Refs. [3 [lOl El 
HH US UHl HTjUH] and references therein.) Recent stud- 
ies have calculated the energy density [53] [3^ and en- 
ergy flux [ST] [52] associated with a heavy quark moving 
through a strongly coupled SYM plasma. These studies 
have shown that qualitative features of the quark wake, 
such as the formation of a Mach cone for supersonic mo- 
tion, match what is expected from hydrodynamics. How- 
ever the comparison between hydrodynamics and the ex- 
act result for the stress-energy tensor has not yet been 
performed in a quantitative fashion. In particular, the 
range of validly of hydrodynamics has not been addressed 
quantitatively. 

Hydrodynamics is valid only on length (or time) scales 
which are large compared to the mean free path (or 
time) of typical excitations in a fluid. In a weakly cou- 
pled relativistic plasma, the mean free path of quasi- 
particle excitations (quarks or gluons) is parametrically 
large compared to their de Broglie wavelength, imfp ^ 
l/{T\'^ln\-^) > ^do Broglie - l/T 33J. The 't Hooft 
coupling A = g^Nc is the appropriate measure of the in- 
teraction strength. As the size of the coupling increases, 
this separation of scales shrinks, and hydrodynamics 
becomes valid on progressively shorter distance scales. 
When A > 1, a description in terms of weakly interacting 
quasiparticles is no longer valid. For an ultra-relativistic 
non-Abelian plasma in this regime, the minimum length 
scale ^hydro On which hydrodynamics is valid must for- 
mally be of order l/T, as there is no other relevant scale. 
But whether, in practice, one needs fhydro ^ l/T or 
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^hydro ^ 100/T in Order for hydrodynamics to reproduce 
results of the full theory, to a given accuracy, cannot be 
determined without a quantitative comparison. Perform- 
ing such a comparison is one goal of this paper. 

To apply hydrodynamics to our situation of interest, 
in which a moving heavy quark is transferring energy 
and momentum to the plasma, one must first formulate 
appropriate effective sources to be used with the hydro- 
dynamic equations. Just as in other applications of long 
distance effective field theories, this requires matching 
the effective theory to the underlying microscopic theory 
at the required level of precision. We show, on general 
grounds, that the appropriate effective sources for hydro- 
dynamics can be simply expressed in terms of the drag 
force acting on the quark. Using gauge/string duality, 
we then compute the energy density and energy flux (or 
momentum density) associated with a heavy quark mov- 
ing through a strongly coupled N — A SYM plasma, and 
compare the exact results to those obtained with hydro- 
dynamics. Comparisons are made in both momentum 
space and position space. We find that the hydrody- 
namic approximation to the quark wake agrees with the 
exact result remarkably well even at distances less than 
2/T away from the quark. Although our comparison ap- 
plies to an A/" = 4 SYM plasma (in the limit of large A and 
large N^) the remarkably good performance of hydrody- 
namics in this context can only bolster the hope that it is 
sufficient to use hydrodynamics to model the transport 
of energy and momentum lost by high energy particles 
traversing a real quark-gluon plasma |34L 135] (despite the 
fact that gradients in the hydrodynamic variables can be 
disconcertingly large) . 

The outline of the remainder of this paper is as fol- 
lows. In Section |IT] we introduce various definitions and 
conventions that we will use throughout the analysis. In 



Section III we discuss the hydrodynamic description of 
perturbations in the SYM stress-energy tensor. This in- 
cludes the formulation of effective sources for the hydro- 
dynamic equations and the separation of the energy flux 
into contributions from sound and diffusion modes. In 
Section |IV| we turn to the dual gravitational formula- 
tion. Gauge/string duality maps the problem of com- 
puting the perturbation in the stress-energy due to the 
quark, /^T^^{x), into the problem of computing the per- 
turbation to the geometry caused by an open string (dual 
to the heavy quark) moving through a five dimensional 
anti-de Sitter/Schwarzschild spacetime. We introduce a 
convenient set of gauge invariant variables which encode 
the metric perturbation, show how their use allows one 
to decouple completely the linearized Einstein equations, 
and then, by analyzing the on-shell gravitational action, 
show how to reconstruct the expectation value of the 
SYM stress-energy tensor using our chosen gauge invari- 
ant variables. We also describe our technique for solving, 
numerically, the ordinary differential equations satisfied 
by the gauge invariant variables, which is based on the 
construction of appropriate Green's functions from nu- 
merically computed homogeneous solutions. We present. 



in SectionjV] the results of the computation of the pertur- 
bation in the energy density and energy flux, in position 
space, due to a moving quark, at several different veloci- 
ties. We also examine the small momentum asymptotics 
and compare results, in several ways, with hydrodynam- 
ics. Section [VT| discusses the interpretation of our results. 
We conclude in Section [Vll| Two appendices contain de- 
tails of the analysis of the boundary action, and the ex- 
traction of small momentum asymptotics. 

We have endeavored to make the presentation rela- 
tively self-contained. Readers interested in the results 
but not the details of the gravitational calculation should 
feel free to skip over Section |IV| 



II. DEFINITIONS AND CONVENTIONS 

We use the Minkowski space metric 77^^ = 
diag(— 1, 1, 1, 1). Five dimensional AdS coordinates will 
be denoted by Xm while four dimensional Minkowski co- 
ordinates will be denoted by x^. Upper case Latin in- 
dices M, N,P, . . . run over 5d AdS coordinates, while 
Greek indices run over the 4d Minkowski space coordi- 
nates. We choose coordinates such that the metric of the 
AdS-Schwarzschild (AdS-BH) geometry is 



-f{u) + dx^ 



du" 
fin) 



(2.1) 



where /(u) = 1 — {u/uh)'^ and L is the AdS curvature 
radius. The coordinate u is an inverse radial coordinate; 
the boundary of the AdS-BH spacetime is at u = and 
the event horizon is located at m = Uh, with T = {■KUh)~^ 
the temperature of the SYM plasma. 

When introducing a Fourier transform over the Ad 
Minkowski space coordinates, we will often decompose 
vectors and tensors in terms of an orthonormal basis of 
polarization vectors {q, Ci, £2}. The quark's velocity v 
defines a preferred direction so a natural choice of polar- 
ization vectors is 



ei 



qx {v X q), €2 



vxq, 



(2.2) 



where q = \q\ and — \q — {ii ■ q) v\ is the magnitude 
of the component of q orthogonal to v. Lower case Latin 
indices a,b, - ■ ■ = 1,2 will be used to refer to the trans- 
verse spatial components of vectors and tensors. We will 
decompose the Fourier transform A^^ {uj, q) of any given 
vector field as follows. 



A^ = A^ 

A' = gM' + e\ A" , 



(2.3a) 
(2.3b) 
(2.3c) 



with a sum over repeated indices implied. The com- 
ponents = {A° , A" , A'' , A^} will be referred to as 
the components of the vector field A in the polarization 
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frame. Higher rank tensors will be represented by polar- 
ization frame components in the analogous fashion. 
It will also prove convenient to use the notation 



Z2 = Zf ea^eb, 



(2.4) 



for quantities which transform with helicity one or two 
under S0(2) rotations about the q axis. 

In four dimensions, any symmetric tensor field T''" 
which satisfies a conservation equation, 



and the trace condition. 



(2.5) 



(2.6) 



for some given and /3, contains only five independent 
degrees of freedom. A convenient representation of these 
independent degrees of freedom is provided by the follow- 
ing helicity zero, one, and two components of the Fourier 
transformed tensor: 



-7- _ -7-00 



r 



(2.7a) 
(2.7b) 
(2.7c) 



where T'^'^ are the components of T'^" in the polarization 
frame. The reconstruction of the original tensor compo- 
nents T'^'^ is given by 

T"° = To , (2.8a) 
= Ti"^ + f , (2.8b) 

+ 1(6^^ -q^q^){P + T^°)^ (2.8c) 

where v^-^u^^ = v^u^ + v^u^ and 

EE {ivTo ~ iV°) /q , (2.9a) 

= {ioT^" - iV") /q , (2.9b) 
T'f' = {ujT°'^ - iV«) /q 

= [lu^To ~ i{iuV° + qV)] /q^ , (2.9c) 

with the components of in the polarization frame. 
We will refer to the quantities {7^} as helicity variables 
and the representation of T^'^ in terms of Tg as the helic- 
ity decomposition of T^" . We emphasize that the helicity 
decomposition is complete only when the helicity vari- 
ables Tg and both (3 and are known. 



III. HYDRODYNAMIC DESCRIPTION 

We consider a fundamental representation quark of 
mass M moving with constant velocity v through an 
A/" = 4 SYM plasma at temperature T. We assume that 



both the quark mass M, and its kinetic energy, are large 
compared to T.^ 

We assume that the quark has been moving at the 
velocity v for an arbitrarily long time At — > 00. Due to 
its interaction with the SYM gauge (and scalar) fields, 
the quark will perturb the surrounding plasma and the 
plasma will exert a friction force, or drag, on the quark. 
The drag force on the quark is minus the rate at which 
momentum is transferred from the quark to the plasma. 
In the absence of any external force, the drag exerted by 
the plasma would cause the quark to lose momentum and 
slow down. To maintain a constant velocity, energy and 
momentum must be supplied to the quark via an external 
force which exactly counterbalances the plasma drag (at 
a given terminal velocity) . This is naturally accomplished 
by turning on a constant background U{1) electric field 
which couples to the fundamental representation quark 
but not to the adjoint representation SYM degrees of 
freedom. 

The microscopic energy- momentum conservation equa- 
tion takes the form 



(3.1) 



where T^'^ is the stress-energy tensor for the system (not 
including the background U{1) electric field), and is 
the external force density, or minus the drag force density, 
acting on the quark. In the limit of large quark mass, 
the quark can be arbitrarily well-localized and the force 
density will have point support. In this regime we may 
write 



F'^it,x) = r S%x-vt), 



(3.2) 



where = dp'^^^^^/dt is the external force (or minus 
the drag force) acting on the quark, and coordinates are 
chosen so that the quark is at the origin at time t — 0.^ 
For strongly coupled SYM, the magnitude of the drag 
force has been evaluated (using gauge/string duality) and 
one finds [31 0] 
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: , (3.3a) 
(3.3b) 

For distances d = \x—vt\ of order l/T or less, gradients 
in the stress-energy tensor are large, non-hydrodynamic 



^ Having kinetic energy large compared to T implies that Brow- 
nian motion of the quark, induced by thermal fluctuations in 
the plasma, may be neglected Our analysis in the grav- 

itational dual will require a stronger condition on the quark 
mass, M 3> V^T, where A is the (large) 't Hooft coupling. 
Strongly-coupled SYM with massive fundamental hypermulti- 
plets has deeply bound mesons whose masses scale as the quark 
mass M divided by \/A, so the condition M 2> T implies that 
both mesons and quarks are heavy compared to T. 

^ Note that ff' is the four-momentum transfer per unit coordinate 
time, not proper time; the covariant four-force dp^^^^^/dr equals 
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degrees of freedom are important, and the complete mi- 
croscopic theory is needed to compute the transport of 
energy and momentum. But as the disturbance in the 
plasma, induced by the quark, propagates out to larger 
distances, dissipative effects will decrease the size of gra- 
dients. At sufficiently large distances from the quark, 
the transport of energy and momentum can be described 
by neutral fluid hydrodynamics. This is the appropriate 
effective theory for the long wavelength, slowly relaxing 
degrees of freedom in a non-Abelian plasma.^ The rel- 
evant hydrodynamic variables are the locally conserved 
energy-momentum densities T'^^(x). In particular, as we 
next discuss, long wavelength perturbations in the spa- 
tial stress tensor can be expressed entirely in terms of 
these conserved densities. 

Instead of working directly with the conserved densi- 
ties r°^(a;), it is conventional, and convenient, to intro- 
duce the proper energy density e{x) and the fluid four- 
velocity uf^{x). The fluid four- velocity is deflned as the 
velocity of a local reference frame in which the spatial 
momentum density vanishes, and the proper energy den- 
sity is the energy density in this local fluid rest frame. 

We will denote with a bar the components of the stress- 
energy tensor in the local fluid rest frame (at the loca- 
tion a;). If the fluid were in perfect equilibrium, with a 
globally-defined rest frame, then in that rest frame the 
stress-energy tensor would have the form. 



■ — £ ) 
=0 

eq — " ' 



(3.4a) 
(3.4b) 
(3.4c) 



with the pressure p and energy density e related by the 
equilibrium equation of state of the fluid, p = p{e). When 
the fluid is not in perfect equilibrium, the definition of 
the local fiuid rest frame allows one to write 



T°\x) = e{x) , 
T°'(ar)=0, 

T'^{x) = p{x) 6ij + Tij{x) , 



(3.5a) 
(3.5b) 
(3.5c) 



where p{x) = p{e{x)) is the equifibrium value of the pres- 
sure which corresponds (via the equation of state) to the 
local energy density e(x), and all non-equilibrium effects 
are contained in the dissipative contribution to the stress 



tensor, r, 



For sufficiently long wavelength variations in the en- 
ergy and momentum density (or proper energy density 
and fiuid velocity), the dissipative contribution to the 



stress may be expanded in terms of spatial gradients of 
the hydrodynamic variables. This, by definition, is the 
hydrodynamic regime. 

It is straightforward to construct the gradient expan- 
sion of Tij. However this exercise is simplified when (i) 
the underlying theory is conformal, and (ii) the variations 
in the energy and momentum density are parametrically 
small compared to their equilibrium value. In conformal 
theories, such as A/" = 4 SYM, the stress-energy tensor 
is traceless. Consequently, e = 3p and tu = 0. If the 
variations in energy and momentum density are small, 
then one may also expand in powers of departures from 
equilibrium. This is the case in our application involv- 
ing a single fundamental representation quark interacting 
with a large SYM plasma. In the large limit, the 
equilibrium stress-energy tensor of the plasma is 0{N^) 
while the perturbations in T'"' due to the addition of a 
single quark are 0(7V°). Consequently, both the frac- 
tional change in the energy density and the fiuid velocity 
induced by the moving quark are of order 1 /N^ . It fol- 
lows that the equations determining perturbations in the 
stress-energy tensor must be linear in the large limit. 
More specifically, nonlinear terms in the hydrodynamic 
equations of motion are suppressed by powers of N^- For 
later convenience, let 



and 



£{x) = e{x)- Ce 



r{x)^fj{x) = l£{x). 



(3.6a) 



(3.6b) 



denote, respectively, the deviation of the energy density 
from its equilibrium value and the associated deviation 
in the pressure. Linearity of the hydrodynamic equations 
of motion implies that only terms linear in u and £ (and 
their derivatives) will be needed in the gradient expan- 
sion of the stress tensor. With this and the vanishing 
trace condition in mind, the gradient expansion of the 
dissipative stress Tij has the form 



V (V^Mj 



\7jUi 



|<5y V • u) 



+ e{ViVj-lSijV^)£ + -- 



(3.7) 



The coefficient r] is the shear viscosity, while the coeffi- 
cient O characterizes the first correction to viscous hy- 
drodynamics. 

The stress-energy tensor in the local fiuid rest frame 
is related to the stress-energy tensor in the lab frame by 
the local boost (in block form) 



A{x) = 



u{x) 



u{x)^ 



+ Oiu{xf). (3.8) 



Hence, the stress-energy tensor in the lab frame is 



^ In Abelian plasmas, the appropriate effective theory is magneto- 
hydrodynamics as the magnetic field can have an arbitrarily long 
correlation length. But in non-Abelian plasmas, even static mag- 
netic fluctuations of the gauge field develop develop a finite cor- 
relation length. 



hydro V-*-^ — -'eq hydro ^^"^^ ' 



where 



Ttq = diag(eeq,Peq,Peq,Peq) , 



(3.9) 



(3.10) 
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and 



hydro 



+ e (v,v, 



f <5,, V • «) 



(3.11a) 
(3.11b) 

(3.11c) 



(neglecting nonlinear corrections subleading in 1/N^). 
Further derivative corrections to the above formulas only 
enter in Eq. (3.11c I, which is the constitutive relation 



expressing the spatial stress tensor in terms of the hy- 
drodynamic variables. 



A. Equations of Motion and Effective Source 

The equations of motion for the hydrodynamic fluctu- 
ations follow from appl ying the energy-momentum con- 
hydro- However, the hydro- 
;ion of the stress tensor is not valid 



servation equation (|3.1|) to T^"^ 
dynamic decomposi 
in the near zone close to where the source F'^ is non- 
vanishing. In this region, gradients become large and 
truncating the gradient expansion in Eq. (3.7 1 is unjus- 



tified, as is the neglect of all non-hydrodynamic degrees 
of freedom. However, in the spirit of effective field the- 
ory, one may formulate effective sources for the hydrody- 
namic equations, which encapsulate the dynamics in the 
near zone. Doing so, we write the energy conservation 
relation for the hydrodynamic fluctuations as 



^M-^hydro(-^) 



(3.12) 



where J'^ is the effective source. Because hydrodynamics 
is a long distance effective theory, the effective source 
can be regarded as having point support at the location 
of the quark. That is, may be expanded in terms of 
delta functions and their derivatives.^ To leading order 
in gradients we have 



j^o) 6\x - vt) 



(3.13) 



The vector can in principle depend on time. However, 
during the (assumed large) interval of time in which the 
quark is moving at constant velocity, j^g^ will be time 
independent. 

The effective source J^(a;) can be determined by cal- 
culating the long wavelength limit of the stress-energy 
tensor via the full quantum theory, and then matching 



* In other words, the spatial Fourier transform of J" may be ex- 
panded in a Taylor series in the wavevector q. This is completely 
analogous to the situation in electromagnetism where, given lo- 
calized charge and current densities, one may similarly represent 
the charge and current densities as a sum of delta functions and 
their derivatives, and this directly leads to the multipole expan- 
sion for the electric and magnetic fields in the far zone. 



its form to hydrodynamics. We do this in Section [V] and 
find, in particular, that the coefficient of the leading 



term in the expansion (3.131 is simply equal the rate at 



which the quark transfers four-momentum to the plasma. 



so that coincides with the external force (|3.3|). This 



value for the leading term in the effective source can also 
be deduced with the following simple argument, which 
nicely highlights the fact that this relation between the ef- 
fective source and the microscopic drag force emerges 
even though hydrodynamics is not applicable in a near 
zone surrounding the quark. 

Assume that the electric field pulling the quark is 
turned on at time ti, and turned off at some later time 
tf. The total four- momentum transfered from the quark 
to the plasma is 



(3.14) 



Ap'^= / dtj^{v{t)) 



where f^{v) is the external force (3.3 1 (or minus the drag 
force) for velocity v. For times t much later than tj, all 
of the four- momentum transfered to the plasma will have 
been transported to distances far from the quark. (As 
discussed below, since the relevant dynamics is diffusive, 
the characteristic distance d is proportional to y/t~tf 
and diverges for large t.) The resulting four-momentum 
stored in the far zone is computable using hydrodynamics 
and, for late times tf, this far-zone four-momentum 



must match the four- momentum (3.14) lost by the quark 



In the limit in which the quark has been moving at 
constant velocity v for an arbitrarily long period of time, 
the total four-momentum transfered to the plasma (ne- 
glecting irrelevant endpoint corrections sensitive to how 
the electric field is ramped up and down) is just 



(3.15) 



In the hydrodynamic description, the total momentum 
transfered to the far zone at time tf is given by 

ApUo(*) = / d'x AT°^^,,,(t, X) . (3.16) 
Rewriting this as 



fjt'^Jd'xAT^^,Jt\x), (3.17) 



and then using the effective energy-momentum conserva- 
tion relation (3.121, yields 



AKyd..o(i) = 



dt' 



d^a 



^« hydro ^ 



(3.18) 



The first term in the integral can be converted to a sur- 
face integral over the sphere at spatial infinity. This sur- 
face term vanishes by causality — if the quark has only 
been moving for a finite period of time, then perturba- 
tions in the stress-energy tensor must vanish at spatial 
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infinity. In the second term, the spatial integral of 
just gives the leading term in the expansion ( |3.13 ) of the 



effective source, / d'^a; — j^^ and, as noted above, j 



'(0) 



will be time independent during the long interval during 
which the quark moves with constant velocity. Conse- 
quently, for t^tf,we have 



(0) 



(3.19) 



Demanding that this reproduce the microscopic result 
(3.151 shows that 



1(0) 



(3.20) 



as asserted above. 

Higher order gradient corrections to , which we eval- 
uate in Section |V] lead to corrections to the stress-energy 
tensor which are suppressed by additional inverse powers 
of I a; I in the far zone. We note however, that in the large 
Nr hmit 



J°(i, x) = F°{t, x) = f -v 5^{x - vt) , 



(3.21) 



with no additional derivative corrections. This follows 
from the fact that the linearized hydrodynamic expres- 
sions for the conserved densities AT^y^j.^, as given in 
Eqs. ( 3.11a[ ) and ( |3.11b I, do not receive any derivative 
corrections (and all non-linear terms are suppressed by 
additional factors of 1/Nc). The form of the gradient 
corrections to J are also constrained. These constraints 
are conveniently formulated in Fourier space. If the du- 
ration of time At in which the quark has been moving at 
constant velocity v is sent to infinity, then J{lu, q) will 
be proportional to 2n5{u! — v ■ q). Furthermore, J((x', q) 
can only depend on the vectors f , q, which implies that 
J{uj,q) must lie in the plane spanned by these vectors.^ 
We therefore can write 

J{uj,q) = [v (j)^{uj , q^) +iq(t>q{uj,q^)\ 2tt5{u:-v ■ q) . 

(3.22) 

The functions (j>v{uj, q^) and 4>q{uj, q^) must be analytic in 
both Lo and q^ , in a neighborhood of = q = 0, so that 
J{t, x) has an expansion in terms of derivatives of delta 
functions with point support at the location of the quark. 
The condition (|3.20| implies that 0.u(O, Q) = f -v/v"^. 



B. Sound and Diffusion Modes 

To compare with the gravitational results presented in 
Sect ion [rv] it will be advantageous to perform a spacetime 



^ In principle one can consider an additional component of J pro- 
portional to V X q. Since this is a pseudovector, its coefHcient 
function must be a pseudoscalar. But there is no pseudoscalar 
that can be constructed out of v and q alone, so this term is not 
allowed. 



Fourier transform and express the stress-energy tensor 
in terms of the helicity decomposition discussed in Sec- 
tion [nj In particular, AT^^^^^ can be reconstructed from 

the helicity variables A7^^^'^'°, ATf'-*"^"^" and AT2^'^'°, as 



summarized in Eqs. (2.8 1-(2.9 1 with V^' = and /3 = 



The ansatz given in Eq. (3.111 and the effective conser 



vation relation (3.121 determine the functional form of 
the helicity variables. Let Jl and Jt be the transverse 
and longitudinal components of J. Ignoring the second 
order Q term and higher order derivative corrections in 
Eq. (3.11c), one finds that the helicity variables are given 

by 



ATn 



hydro 



p 



^^hydro ^ 



Dq^ ' 



^^thydro 



0, 



where 



p^iq 



J + iuj J° -jq'^J^ . 



(3.23a) 

(3.23b) 
(3.23c) 

(3.24) 



Here Cg — de/dp is the speed of sound, 7 = 4?7/3(e + p) is 
the sound attenuation constant, and D = ri/(e + p) is the 
transverse momentum diffusion constant. For strongly 
coupled SYM IMj 



Cs = 1/3, 

7 = i/(3^r) 

D = l/(47rr) 



(3.25a) 
(3.25b) 
(3.25c) 



The lack of any helicity two component, Eq. (3.23c I , re 



mains true (in the large A^c limit) even when higher order 
gradient corrections are included. This is a consequence 
of the fact that ATj^^^^.^ is linear in the perturbations 

£ and u*. With this constraint, regardless of the num- 
ber of spatial derivatives, one cannot construct a nonzero 
transverse traceless component of the spatial stress. 

It is straightforward to reconstruct the perturbation in 
the energy density £ and energy flux Si = ATj^y^j.^ from 
the helicity variables and the effective energy-momentum 
conservation relation (3.12). The energy density and en- 



ergy flux, in spacetime, are given by 



£{t,x) 
SL{t,x) 
ST{t,x) 



dbj (fiq 
27r (27r)-^ uj"^ ^ c^q^ -\- i"f (pu 



P 



_ I" duj d^q cl iqj" + iujJl iQ., 
27T (27r)'^ — q'^ + ij q^u) 

du) d^q — Jt jQ.j, 

2^(27r)3 iuj-Dq^^ 



(3.26) 

■X 

(3.27) 
(3.28) 



with S{t,x) = Si,{t,x) + ST{t,x) and the four-vector 
Q = (cij, q). These expressions show that the energy den- 
sity and the longitudinal component S'l of the energy 
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flux satisfy a diffusive wave equation, while the trans- 
verse component St of the energy flux satisfles a diffusion 
equation. The energy density obeys the inhomogeneous 
wave equation 



(3.29) 



which describes damped sound waves, traveling at speed 
Cs — l/-\/3, with a source p given by Eq. (3.241. 



It is instructive to separate the energy flux into sound 
and diffusion modes. Since the longitudinal part of the 
flux 5'l satisfles a wave equation while the transverse 
part St satisfles a diffusion equation, one would natu- 
rally expect that should be identifled with the energy 
flux carried by sound waves, while St characterizes dif- 
fusive energy flux. But, in a steady-state situation where 
the quark has been moving for an arbitrarily long time, 
thin gs are not so simple. The effective source J{uj,q), 
Eq. (3.221, equals 5{oj — v ■ q), times a function which is 
regular as uj and q — )■ 0. But its longitudinal projection 
has a directional singularity at 9 = 0, 



Jl = 



= q 



2tt5{lo—v ■ q) , 



(3.30) 

with • • • denoting terms regular as g ^ 0. Therefore, the 
Fourier transform of is not localized at the quark, but 
rather falls off in space like the inverse cube of the dis- 
tance away from the quark. As a res ult, the longitudinal 
energy flux Si^, deflned by Eq. (3.271, contains the term^ 



nlocal 



dbJ dPq 
2^(27r)3 ^ 



iq 



q^) 2-k5{lo—v ■ q) e 



iQ-x 



(3.31) 



This gives a Coulomb-field-likc flux at large distance, 



*^nonlocal(^7 ^ ^ 



1 



AttvIx — vt\ I 



f v 



(3.32) 



This is co-moving with the quark and does not propagate 
like a wave. It makes little sense to regard this as a contri- 
bution to the energy flux carried by sound waves. (Note, 
in particular, that ■S'noniocai is non- vanishing in front of 
the quark even when the quark is moving supcrsonically.) 



^ To see this, use the frequency delta function to rewrite the factor 
of UJ multiplyi ng J\ in Eq. (|3.27|l times the {v ■ q) factor in the 



c^q + i'^q^w) + q (c^ — i'yuj) 



singular term ( |3.30| 
The first term cancels the denominator of Eq. ( |3.27| , producing 
the above result for Snoniocal- The second term is no longer sin- 
gular when divided by and is properly viewed as a contribution 
to the energy flux carried by sound. In the resulting expression 



( 3.31 1, if one writes 4>viuJ,q'^) = 4>v(u}, 0)— \4>v{ijJ, q^) — <f>v{uJ, 0)] 



then, strictly speaking, it is only the first </)„(a;,0) term which 
gives a non-local contribution. But including the additional 
\ij>v{uj, q^) — <f>v(uj, 0)] term in the definition of S'noniocai leads to 
simpler expressions for the effective sound and diffusion sources 
discussed below. 



Since the transverse part of the effective source, Jt 
equals J (which is regular as q — > 0) minus Jl, a com- 
pletely parallel argument shows that the transverse flux 
St, as deflned in Eq. (3.28 1, contains a non-diffusing con- 



tribution which is precisely —S'noniocai- 

The presence of these equal and opposite Coulomb- 
like contributions naturally suggests that the energy flux 
which is properly associated with sound and diffusion 
can be obtained by suitably adding and subtracting this 
term. To this end, we write 



*sound 



'diffusion i 



with 



*^sound — "^L '^nonlocal i 
^'diffusion = St + S'noniocai ■ 



(3.33) 



(3.34a) 
(3.34b) 



Using Eq. (3.271, it is straightforward to see that S'sound 



satisfles the damped wave equation 

{~di + cl V + j\/^do) Ssound = Jsound (3.35) 

with a source 

Jsound{t,x) = J ^-^^^Q {2n)5{uj - V ■ q) e''^'^ 

X [c^f + {cl-i-fuj)(j)^{uj,q^) + iu;(l)g{u;,q^)] . 

(3.36) 

Unlike Ji^{u},q), the integrand dcflning the modifled 
source Jsound is now regular as q ^ 0. Hence Jsound 
is localized at the quark (it may be expanded in terms 
of delta functions and their derivatives) . Similarly, using 
Eq. (3.281, one flnds that ^diffusion satisfles the inhomo- 



geneous diffusion equation 



{do - DV^) Sdi 



ffusion 



= J: 



diffusion : 



(3.37) 



with the source 



Jdiffusion = / !!^77r-|r (27r)5(w-v • q) e^^'' 
X {v + iqD)(t)v{uJ,q^) , 



(3.38) 



which is also localized at the quark. Because the re- 
deflned sources Jsound and Jdiffusion are localized, the 
corresponding energy fluxes, S'sound and ^diffusion, may 
be regarded as a sensible decomposition of the flux into 
contributions from propagating and diffusing degrees of 
freedom.^ 



Since V ■ S'noniocai is localized at the quark, the fluxes S^ound 
and Sdiffusion do provide a decomposition of the energy flux into 
pieces which are, respectively, longitudinal and transverse every- 
where in space excluding the position of the quark, up to terms 
which fall faster than any inverse power of distance from the 
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At leading order in the gradient expansion of J, the 
sound and diffusion sources are given by 



Jsound = c2 (f + VF" , 

*^diffusion — ^ i 



(3.39a) 
(3.39b) 



while the leading behavior of the source (3.241 for the 
energy density is 



p^V -F- OqF^ 



(3.40) 



Using these sources and the corresponding wave and dif- 
fusion equations (3.29), (3.351, and ( 3.37[ ), the long range 
behavior of the sound and diffusion modes is completely 
specified. We note that in the large limit the magni- 
tude of the drag force appears as an overall normalization 
of the perturbation in the stress-energy tensor. 



IV. GRAVITATIONAL DESCRIPTION 

According to gauge/string duality. A/" = 4 SYM (on M^) 
is equivalent to type IIB string theory on AdSc, x 
[18lll9j . Turning on a non-zero temperature corresponds, 
in the dual description, to introducing a black brane (a 
black hole with a flat horizon) into the ^^5*5 space, lead- 
ing to t he A dS-Schwarzschild geometry described by the 
metric (lO.^ 



The addition of massive fundamental representation 
fields (an TV =2 hypermultiplet) to the A/" = 4 SYM the- 
ory corresponds, in the dual gravitational description, to 
the addition of a D7 brane wrapping an of the internal 
and covering the five dimensional asymptotically AdS 
space from the boundary down to a minimal radial posi- 
tion Um, which is inversely related to the hypermultiplet 
mass M for large mass [4]. 

A single quark moving through the Af—4 SYM plasma 
corresponds, under gauge/string duality, to an open 
string which runs from the D7 brane down to the black 
hole horizon [H |37] . The presence of the string perturbs 
the geometry via Einstein's equations and the behavior of 



quark. Because hydrodynamics is only valid sufficiently far from 
the quark, one may regard Sgound ^^'^ ^diffusion only being 
defined in M'^\B, where B is a ball surrounding the quark. In 
such a non-contractible region, the decomposition of a vector 
field into transverse and longitudinal components is not unique. 
The above decomposition is a choice for which the resulting effec- 
tive sources Jsound a-iid ^diffusion are localized at the quark (i.e., 
have Fourier transforms analytic in q). With that additional 
condition, this decomposition is effectively unique — alternative 
choices [such as the replacement of (/)t,(aj,q^) by (/)„(aj,0) in the 
definition ( |3.31| of Snonlocal] merely correspond to adding to 
iSsound (and subtracting from iSdiffusion) a- localized contribution 
which falls exponentially with distance from the quark. 
® Turning on a temperature does not deform the internal S^. Only 
the AdS-Schwarzschild spacetime will be relevant for the our pur- 
poses; the five-sphere will play no role and may be ignored. 



the metric perturbation near the AdS boundary encodes 
the change in the SYM stress-energy tensor.^ 

In the Nc ^ oo limit, the 5d gravitational constant 
becomes parametrically small and consequently the pres- 
ence of the string acts as a small perturbation on the 
AdS-BH geometry. To obtain leading order results in 



N[/Nc <C 1, we write the full metric as G 



MN 



G 



IMN, 



where G^-^jj^ is the metric of the AdS-BH geom- 
etry given in Eq. ( |2.1[ ), and then linearize the resulting 
Einstein equations in the perturbation Hmn- 

According to the AdS/CFT correspondence, the on- 
shell gravitational action Sq is the generating functional 
for the boundary stress-energy tensor [38] [39l |40]. The 
metric Gmn induces a metric g^j^i, on the boundary of the 
AdS-BH geometry [ID]. The boundary metric is related 
to the boundary value of G^^ by a scaling function of 
the AdS radial coordinate. More specifically, because the 
metric Gmn has a second order pole at the boundary, the 
boundary metric, which must be regular at the boundary, 
may be defined as 



(0) 
MN 



= — r 



(4.1) 



where the infinitesimal e will be sent to zero after the 
required boundary terms needed to properly define the 
gravitational action are added. The expectation value of 
the boundary stress tensor is then given by [HJ HO] 



T^^ix) = lim ^ ^ 



(4.2) 



with g denoting the determinant of g^jy. Defining for 
later convenience 



(4.3) 



the boundary stress-energy tensor may also be expressed 
as 



T^"'{x) = 21im 



5Sg 



SH^^{x,e) ' 



(4.4) 



(Since the boun dary is flat Minkowski space, the factor 
of ^/—g in Eq. (4.2 1 approaches one as e — > and may 
be ignored.) 

The action for the entire gravitational system is given 

by 



EH 



GH 



s- 



DBI 



5, 



NG 



CT ■ 



(4.5) 



® The D7 brane also perturbs the geometry. This is a small correc- 
tion of order 1/Nc provided the number of flavors N{ (which is 
the same as the number of D7 branes) is held fixed as A'c — > oo. 
This correction encodes the 0{Nc) contribution of fundamen- 
tal representation fields to the equilibrium pressure and energy 
density of the SYM plasma. We ignore this correction (which is 
subleading relative to the 0{N^) contribution of adjoint repre- 
sentation fields) as our goal is the non-equilibrium contribution 
due to a moving quark. 
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The first term is the Einstein-Hilbert action, 



^1 



EH 



1 

2nl 



d^x (i?+2A) , 



(4.6) 



where Kg 
G is the determinant of G^i, and A 



Att'^L^/N^ is the 5d gravitational constant, 

is the cosmo- 

logical constant. The next term is the Gibbons-Hawking 
action. 



-7 2if . 



(4.7) 



with 7pj, the induced metric on the slice u = e (and 7 
its determinant), and K the trace of the extrinsic curva- 
ture of the slice. This term is needed to obtain an action 
which makes the Dirichlet problem well defined (i.e., an 
action depending only on first derivatives of the metric) 
[41j . The remaining terms are the Dirac-Born-Infeld ac- 
tion 5'dbi for the D7 brane, the Nambu-Goto action S'ng 
for the string, and a counter-term action ^ct defined 
on the boundary which is required for holographic renor- 
malization |40j . This term cancels poles in 1/e which 
are generated by the bare gravitational action; the form 
of the counter-terms is constrained by the requirement 
of diffeomorphism invariance of the gravitational theory. 
In Eq. (4.5 1 all boundary integrals are evaluated on the 
u — e slice, with e ^ after all terms are combined. 



A. Heavy Quark Effective Theory 

In the limit that the quark mass becomes arbitrarily 
large, its presence in the plasma merely serves as a exter- 
nal source for the SYM fields. This can be made explicit 
by constructing a heavy quark effective theory (HQET) 
from the field theory Lagrangian (as given explicitly in 
Ref. 113). As we now outline, this procedure has a nat- 
ural counterpart in the gravitational dual. 

The interaction between the string and the D7 brane 
is governed by the DBI action and the Nambu-Goto ac- 
tion. Fig. [1] shows a cartoon of the string plus D-brane 
system. As the string moves, a fiux of energy and mo- 
mentum flows down the string to the black hole horizon 
[3J H]. This flux of energy and momentum is responsi- 
ble for the drag force on the quark in the dual boundary 
theory and is supplied by a U{1) gauge fleld living on 
the D7 brane — the strength of which is tuned to match 
the drag force acting on the quark and thereby maintain 
a constant velocity. The presence of the string pulls on 
the D7 brane, deforming it. However, in the large mass 
limit the trailing string will only deform the D7 brane 
over length scales of order 1/M. As M — > cx), the string 
endpoint approaches the boundary and the size of the 
region in which the D7 brane is significantly deformed 
shrinks to zero.^*^ 




Black Hole 



FIG. 1: A cartoon of the string plus D7 brane system in 
the large mass limit. The D7 brane covers the asymptotically 
AdS space down to a minimal radial position, away from the 
string, denoted Um- The trailing string is moving to the right 
at constant velocity v. An energy flux flows down the string 
toward the black hole horizon, which is located at radial coor- 
dinate Uh- This energy is supplied by a constant U{1) electric 
field living on the D7 brane. The D7 brane is deformed in 
the neighborhood of the endpoint of the string over a length 
scale of order 1/M. 



Instead of solving the string plus D7 brane system in 
its entirety, just as in HQET one may consider the large 
mass limit and construct an effective gravitational theory 
for the perturbations in the SYM stress-energy tensor. 
In the effective theory we send the quark mass M 
00 before performing holographic renormalization of the 
gravitational action. In doing so, we take u„j — > so the 
D-brane no longer extends inward from the boundary. 
Furthermore, the deformation of the D-brane due to the 
string shrinks to a vanishingly small region about the 
string endpoint (which approaches the boundary) . 

This suggests that the dynamics of the D7 brane should 
be irrelevant in the heavy quark limit. This is not quite 
true, however. Because general relativity is a gauge the- 
ory, the construction of a heavy quark effective theory is 
constrained by the requirement that the effective theory 
be invariant under infinitesimal diffeomorphisms 



Xm Xm + £,m{X) , 



(4.8) 



This is a bit oversimplified. For a finite mass M, the back- 



where S^M is an arbitrary infinitesimal vector field. Be- 
cause an energy-momentum flux flows from the boundary 



reaction of the string on the D7 brane is expected to cause the 
D-brane embedding to develop a thin narrow tube which reaches 
all the way down to the horizon 1431 144| . In other words, the fun- 
damental string will be "puffed up" into a tiny tube. However, 
the resulting dynamics of a sufficiently thin tube is indistinguish- 
able from that of a fundamental string. In the M — > 00 limit, 
the width of the tube (at any non-zero value of u) will shrink to 
zero and this issue may be ignored. 
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down the trailing string, gauge invariance at the bound- 
ary requires that the U{1) electric field on the D7 brane 
not be neglected. We therefore construct the effective 
theory for the perturbations in the SYM stress-energy 
tensor with the brane effective action 



B. Gauge Invariants 



rroff 



DBI 



Sem , 



(4.9) 



where Sbm is the Maxwell action for the U{1) electro- 
magnetic field which resides on the boundary. 

To complete the definition of the effective theory, we 
must specify the counter-term action. The counter-terms 
must cancel 1/e poles in the gravitational action and 
must be diffeomorphism invariant. The presence of the 
string ending on the boundary induces 1/e poles in the 
gravitational action which have point support at the lo- 
cation of the quark. The remaining 1/e poles in the grav- 
itational action are simply those obtained in pure gravity. 
The appropriate counter-term needed to offset the pure 
gravity divergences is well known 140 L I45j and given by 



ScT 



1 

'24 



(4.10) 



This term is simply proportional to the area of the bound- 
ary. 

Point-like 1/e divergences with support at the loca- 
tion of the quark are fundamentally different from the 
other UV divergences — they are physical as the bound- 
ary stress-energy tensor should contain a 1/e divergence 
at the location of the quark. This simply reflects the 
contribution to the stress-energy tensor of an infinitely 
massive quark. 



(4.11) 



where is the quark's four velocity. The (divergent) 
coefficient M can either be determined by a matching 
to the divergent point-like terms in the boundary action 
or, alternatively, by integrating the energy density of the 
trailing string and isolating the divergent contribution to 
find its bare mass. Taking this approach, we find |1] 



M = 



27re 



(4.12) 



Instead of adding local counter-terms to the gravitational 
action to offset the point-like divergence, we simply define 
the boundary stress-energy tensor via 



T'^'ix) = lim 



6Sr 



^quarkl-^/' 



(4.13) 



with M and e related by (4.12). The resulting stress- 
energy tensor is finite as e — *■ and everywhere traceless, 
including at the position of the quark. 



So the Fourier transform T^"'(a;, q) of the stress-energy tensor is 
traceless for all momenta. 



The continuity equation (3.1 ) and vanishing trace con- 



ditions satisfied by the M — A SYM stress-energy ten- 
sor imply, as discussed in section |IT] that the stress- 
energy tensor contains five independent degrees of free- 
dom which are conveniently isolated by performing a 
spacetime Fourier transform and decomposing the stress- 
energy tensor in terms of helicity variables Tq, Ti and T2. 

The five independent degrees of freedom contained in 
T^i, may be contrasted with the 15 degrees of freedom 
contained in the metric perturbation Hmn- Not all of 
these degrees of freedom are physical, however. The 
linearized gravitational field equations are invariant un- 
der the infinitesimal diffeomorphisms (4.8 1. Under such 



transformations, the metric perturbation transforms as 



l-MN 



Dm iN — Djy ^ 



M 



(4.14) 



where Dm is the covariant derivative with respect to the 
background metric c'f^jq. Physical degrees of freedom 
carried by Hmn must be invariant under the above gauge 
transformations. This limits the number of independent 
physical degrees of freedom carried to the boundary by 
the metric perturbation Hmn to five, matching that of 
the SYM stress tensor [25] . 

The correspondences between the number of indepen- 
dent gauge invariant degrees of freedom contained in the 
metric perturbation Hmn ^^nd the number of independent 
degrees of freedom in the SYM stress tensor suggests that 
the bulk to boundary problem can be formulated directly 
in terms of gauge invariant degrees of freedom. Using 
this approach has several advantages |35]. As we show 
below, the gauge invariants can be chosen to satisfy de- 
coupled equations of motion. Furthermore, as we show 
in Section |IV D| the on-shell gravitational action can be 
expressed in terms of gauge invariant degrees of freedom 
plus computable boundary terms (which are due to the 
non-conservation of the SYM stress tensor at the location 
of the quark). 

Gauge invariants can be constructed out of the Fourier 
mode amplitudes hMN{u;oJ,q) and classified according 
to their helicity under rotations about the q axis [46] . As 
we shall discuss in detail below, the determination of the 
SYM stress-energy tensor requires the construction of a 
scalar gauge invariant Zq, a vector gauge invariant Zi, 
and a tensor gauge invariant Z2. The choice of gauge 
invariants Zg is not unique. However, all gauge invari- 
ants of a common helicity carry the same information 
to the boundary. We are therefore free to chose any 
convenient set of gauge invariants. 



This may be verified by analyzing the asymptotic behavior of 
both the linearized Einstein field equations and the various gauge 
invariants of the same helicity near the boundary. Doing so, one 
finds that the boundary values of gauge invariants (with the same 
helicity) are linearly related. 
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Let H-MN denote the components of Hi\in in the po- where 



larization frame. Using Eq. (4.14) it is easy to see that 



the following five degrees of freedom are gauge invariant 



Zo = U 



00 -T 2ujqHoq 

Z2 = {Hab - ^"Hcc Sab) ^a®h, 



(4.15a) 
(4.15b) 
(4.15c) 



where sums over repeated indices are implied and primes 
denotes differentiation with respect to u}^ 



C. Equations of Motion 

The equations of motion for the invariants Zs follow 
from the linearized Einstein field equations. The full Ein- 
stein field equations are 



R 



MN 



\Gmn [R + 2A) — Kg Imn ■ 



(4.16) 



where Imn is the bd stress-energy tensor of the trailing 
string. Writing 



G 



MN 



^ MN "MAT : 



(4.17) 

ng the 

left hand side of the field equations (|4.16|) in the pertur- 



where cj^^jy is the AdS-BH metric, and expanding the 



bation Hmn, produces the linearized equations, 

— Hmn + "^D^ D(MhN)P — DmDn h+ h-MN 



{D^h - D^'D'^ hpQ -^^h) g'^j^ = 2kI tMN , 



(4.18) 



where h = h^^. 



The trailing string profile is determined by minimizing 
the Nambu-Goto action for a stationary, constant veloc- 
ity profile. The result is 



X{t,u) ^Vt + a;string(M) , 



with 

•^string ('^) 



VUh 



1 / W 

tan~M — 



:ln 



Uh-U 



(4.19) 



(4.20) 



The 5d stress-energy tensor for the trailing string is |26j 
ioo - s (/ + i;2^,4^-4) ^ t55-s(/ + i;2)r2, (4.21a) 
hi = -svi, tij = sv^Vj , (4.21b) 

to5 = -s r'^u^uj;^ , U5 = s V, r^u^u^"^ , (4.21c) 



Using the differential helicity one invariant l |4.15b| instead of an 
equivalent non-differential invariant turns out to be more con- 
venient for technical reasons involving the numerical stability of 
the resulting ordinary differential equation. 



s{u) 



u \/A 



27ri3%/T^ 



5^{x-vt~x^t,;^^{u)) . (4.22) 



We denote the components of t^^ in the polarization 
frame by t^^. 

Because the background geometry is invariant under 
spatial rotations, the differential operators on the LHS 
of Eq. (4.18 1 cannot couple linear combinations of TLmn 



which transform with different helicities under rotations 



about the q axis. Consequently, the field equations (4.18 1 



reduce to three coupled sets of ordinary differential equa- 
tions labeled by their helicity. Although it is not neces- 
sary, working in the gauge in which /15JV = for all N 
simplifies the extraction of decoupled equations for the 
gauge invariants Zg. The scalar set of equations then 
reduces to four second order and three first order equa- 
tions, while the vector set of equations reduces to two 
second order and one first order equation. However, not 
all of these equations are independent — some of the sec- 
ond order equations can be derived from the first order 
equations. One may use this observation to reduce the 
number of required second order equations of motion. 



1. Tensor Mode 

The components of Ti^j^y which transform with helicity 
two under rotations about the q axis are the traceless part 
of TLab- The ab components of the linearized equations 



(4.181 may be written explicitly as 



-J l^ab l^ab H r l^ab 



f 



(4.23) 



Jab 



tab : 



where Q'^ = (uj/f.q) and 

K^ab = Ti-ab + (Woo/ / — Ti-ii) Sab ■ 



(4.24) 



The equation of motion for the helicity 2 gauge invari- 
ant Z2 immediately follows from the traceless part of 



Eq. (4.231, 



^2 Z'2 + B2 Z2 — 82, 



where 



A2 
B2 
S2 



uf - 3/ 
uf ' 

q^f - 



P 



2Ki 



{tab — \ tec Sab) ^a® 



(4.25) 

(4.26a) 
(4.26b) 
(4.26c) 
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Inserting the st ring stress-energy tensor (4.211 into ex- 
pression (4.26c) for the source 5*2 gives 



27rL3 q^fVl' 



{2Tr)d{LU-v ■ q) 



X e ^-J-^"'""* (ei (g) _ £2 «) £2) 



2. Vector Mode 



(4.27) 



The components of 'Hfn, which transform with hehc- 
ity one under rotations about the q axis are Tioa, Wag 
and 7ia5 (which has been set to zero). The equations 
of motion for these quantities contain two second order 
equations and one first order equation. For the deriva- 
tion of the equation of motion for Zi , it is sufficient to use 
one second order equation (which is independent of the 
first order equation) and the first order equation. We use 
the following components of the linearized field equations 



(4.181 



3/ 



H'oa+q^nOa+q^nag=2Kltoa, (4.28) 



/ 



'icjn'oa + ^qfKg)^2nlta5. (4.29) 



Using these equations, we find the equation of motion for 
the helicity one invariant Zi 



where 



Bi 
Si 



Z'; + AiZ[ + Bi Zi^Si, 

uf - 3/ 
uf ' 

3p - u {uq"" + 3/') / H- u'uj'' 



/ 



(4.30) 

(4.31a) 
(4.31b) 
(4.31c) 



Subs titutin g the string stress-energy tensor (4.211 into 
Eq. (4.31c I, gives 



4Vx 



vqj_ 



Si^^ — ^ {2n)Siu;-v ■ q) e-"^-- 



(4.32) 



3. Scalar Mode 

The (non-vanishing) components of Ti^^ which trans- 
form like scalars under rotations about the q axis are 
WoOi T^Oqi T^qq and Tiaa- The equations of motion for 
these quantities contain four second order equations and 
three first order equations. For the derivation of the 
equation of motion for Zq, it is sufficient to use one sec- 
ond order equation (which is independent of the three 



first order equations) and the three first order equations. 
By talcing suitable linear combinations of the components 
of the linearized field equations (4.18 1, we find the follow- 
ing equations of motion: 



nill 

'T-oo 



f^^f nil I / nil 

-^^00 + jn,,- 

3/ " 3 



12 



6/2 



■n 



00 



Aqcu 



Oq 



3/ 



2tn 



(4.33a) 



f 



00 



/ 

2uLoq 



Oq 



Uq Uaa 
:2ufKlt55, (4.33b) 



iuH'og + iqHoQ-iqfHa 



iujU'.i + iqU'Qq - 



2/ 



'Ha — 



^qf^ 
2/ 

iqf 



(4.33c) 



HQq 2Kc; ti 



5 E05 



(4.33d) 



It is straightforward (but tedious) to derive the equation 
of motion of Zq from the above equations. The result has 
the form 



Zq + Aq Z'q + Bq Zi — So . 



(4.34) 



The coefficients and Bq and the source Sq can be 
determined by substituting the definition of Zq into the 
above ansatz and then using the above equations of mo- 
tion (4.33a|-(4.33dl to eliminate all derivatives except 



TYqq. Demanding that the coefficient TYqq vanish then 
determines ^o- Once Aq is determined, it is straightfor- 
ward to read off the coefficient Bq and the source 5*0. 
One finds 



Ao 



uf 2i{q^f-u;^) 
f + q2 („j,_6/) + 6a;2 



/ q^uf'-6f) + 6Lo2 



(4.35a) 
(4.35b) 



The expression for the source Sq in terms of a general 
string stress-energy tensor tpiN is sufficiently lengthy 
that we will not give it here. But substituting the explicit 
form of the string stress-energy tensor (4.211 produces a 
relatively compact result. 



Sq 



kIVX q^ {v^+2) -SuJ^ 



67rL3 q2yX^ 

u[q'^u^ + A&iq^uulu^ - Q{q^-uj'^Yul] 
/(/g2 + 2g2_3L.2)u8 

X {2ti)5{uj-v ■ q) e-'i-"'^":""^ . 



(4.36) 
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D. SYM Stress Tensor 



The variation of the on-shell gravitational action may 
be expressed as 



5Sg = SS] 



hori: 



5Sb 



(4.37) 



where Sb is a surface term at the boundary, and ^horizon 
is a surface term at the horizon. Following Ref. [47] , we 
neglect ^horizon ■ We show in Appendix [A| that 



{2ny 



Ao SZl du + Ai 5Z\ Zi 



A2tr 5ZlduZ2 + \5Hl,%''^^ 



-SS^ 



EM 



0(e) 



ith 



Aq 

Ai 

A2 



1 

1 



(4.38) 

(4.39a) 
(4.39b) 
(4.39c) 



and 



(4.40) 



In the fourth term of (4.381, Tj^^ denotes the stress- 
energy tensor of the equilibrium N — A SYM plasma. 



7^cT= i^cVT'i diag(l,i,i,i). 



(4.41) 



with its components in the polarization frame. In 
the fifth term of (4.381, the symmetric tensor 3'*'' is a 



linear combination of components of the string stress- 
energy tensor evaluated at the boundary. Explicitly, 



This can be understood as follows. As we discuss below, causality 
implies that near the horizon the gauge invariants behave as 
Zs ~ (m—u^ To make the u —> ui,^ limit of this quantity 

(and Shorizon) meaningful, the frequency can be infinitesimally 
analytically continued into the complex plane. Causality dictates 
the required direction of continuation. For retarded boundary 
conditions, one must send ui into the upper half plane, lj — > u+iS. 
With such an infinitesimal continuation understood, the gauge 
invariants Zs v anish at th e horizon. 

See Eqs. | |A27| and |A3i} for the definitions expressing 3^" in 
terms of tjvfiv- 



311 = 
= c 

with 




icu [g2 (5^2+1) - 3c^2 („2^^)] 
iu; [3g2 („2+i) _ („2+5) ^2] 

U2 (l-^;2) (g2_^2) 

qu; ^[3^;2g4+(l„„2)^2g2_3^4 



(4.42a) 



(4.42b) 



quf^ (1—^2) (q'^—u^) 



q'^—uj'^ 



lUJ 

7/2 



3ivq± {(f—uj'^) 
q^uj^ (1— f 2) 



C 



(27r)(5(w-v • q). 



(4.42c) 
(4.42d) 

(4.42e) 



(4.43) 



All other components of 3^'^ vanish. In Eq. (4.42d), no 
sum on a is implied. 

On the boundary, the action of diffeomorphisms on the 
metric perturbation takes the simple form 



(4.44) 



The first four terms in the boundary action (4.38) are 



manifestly invariant. The fifth term however violates dif- 
feomorphism invariance. One may easily verify that 3^'^ 
satisfies 



where 



^0 



(gg ■ v)V\ 



{2Tr)5{uj-v ■ q) , 
{2tt)5{lo-v ■ q) , 
{2tt)5{uj-v ■ q) . 



(4.45) 

(4.46a) 
(4.46b) 
(4.46c) 



These are precisely the components of the external force 
density F'^ in the polarization basis. The violation of dif- 
feomorphism (or gauge) invariance only occurs at the lo- 
cation of the string endpoint. The presence of 3^" in the 
boundary action refiects the fiux of energy-momentum 
from the boundary into the string. This fiux comes from 
the U{1) electric field living on the D7 brane. By defini- 
tion. 



5S 



EM 



(27r) 



4 2 



(4.47) 
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where T^j^ are the components of the electromagnetic 
stress-energy tensor in the polarization frame. By con- 
struction, the strength of the C/(l) electric field is tuned 
to offset the drag force and maintain a constant velocity 
of the quark. It therefore follows that we must have 



EM 



(4.. 



Consequently, the total boundary action is gauge (diffeo- 
morphism) invariant. 



Given the boundary action (4.38 1, one can evaluate the 



boundary stress-energy tensor via the definition (4.131 



As discussed earlier, it is convenient to express the SYM 
stress-energy tensor in terms of its helicity variables 7^ . 
Let r^q^ark j^^^ ^^cq denote the helicity variables of T^^^^^ 
and T^!^ , respectively. Neglecting the contribution from 
the C/(l) electromagnetic field, a short exercise yields 



To ~ lim 
Ti = lim 
T2 = lim 



Zq -r o ■'0 
2A2 duZ2 



n 



uark 



, (4.49a) 
(4.49b) 
(4.49c) 



We now simplify the above expressions for the helicity 
variables Tg. The gauge invariants arc solutions of their 
respective differential equations (4.251, (4.30) and (4.341, 



and vanish at the boundary. This implies that they have 
the power series expansions. 



Za 
Zi 



u'Z^^ 

3*5(3) 



Z2 = u'Z^ 



'u'zt^ 



(4.50a) 
(4.50b) 
(4.50c) 



The leading terms in these expansions are temperature 
independent. Plugging the expansions into the respective 
equations one finds 



Z, 



(3) 



4V\ [q^ -3c^2] (^2_^2) 



Zl 



(2) 



67rL3 



X (27r)5(a; — V ■ q) , 



vqj_ 



z: 



(3) 



? 9 



67rL3 g20-3^ 



X (ei 



{2tt)6{uj—v ■ q)€i 



{2Tr)6{uj-v ■ q) 
•h-h® £2) ■ 



(4.51a) 
(4.51b) 

(4.51c) 



From the definition (4.111 of T^^^^^ and our choice (2.2) 
of polarization vectors, it is easy to see that 



T 



quark 



M 







{2n)5{uj^v-q), 



quark 



/^-quark 



M 



vq±_ 



(2'k)S(lu~v ■ q) ei 

[vqi-) 



2 q 
M ^^.^.^2 



-V 



2g2 
X (ei 



{2tt)5{uj-v ■ q) 
' ei - £2 <8) £2) , 



(4.52a) 
(4.52b) 

(4.52c) 



where M = •\/A/(27re) is the quark mass which is being 
sent to infinity. The helicity components of the equilib- 
rium stress-energy tensor are just 



8 

Putting it all together we have^^ 
4g2i3 



To 



4"' + I? + ro°^ 



3K2(g2_^2)2 
^^(3) 

5 



1 2 - —2- Z2 



(4.53) 

(4.54a) 
(4.54b) 
(4.54c) 



where 



iuj (5w2-f 1) - 3 (w^-1) uj^ 

1*2 (1-1,2) (g2_^2) 



\/1—v^{2'k)S{uj—v ■ q) 

iLo [g2 (5«2+l) - 3 {v^-l) u;2] 

ul (1-^2) (g2_^2)2 



(4.55) 



Note that, after subtracting 7^'^"'*'''*^, all divergent quanti- 
ties have canceled. 

Below, we will focus on the temperature dependent 
perturbation in the stress-energy tensor due to the mov- 
ing quark, 



T = 



(4.56) 



That is, we subtract the zero temperature stress-energy 
of the quark as well as the stress-energy of the equilibrium 
plasma. We also define 



= AT°°, 
A5, = AT°\ 



(4.57a) 
(4.57b) 



which are the temperature dependent perturbations in 
the energy density and energy flux, respectively. 



^® These expressions for the Ts can also be obtained by solving 
the complete set of field equations for the metric perturbation 



Near the boundary, H^^v 



H 



(3) 



(4) 



In 



a gauge in which /ism = 0, the SYM stress-energy tensor is 
given by T^j/ = + Hj;^J [40) . This approach was taken 

in Refs. 1261 1291 1481 . Our treatment above has the virtue of 
being self-contained and very explicit about the treatment of 
contributions from the quark and the external electric field. 
^'^ The zero-temperature stress energy tensor is completely deter- 
mined, up to the overall normalization, by conformal invariance. 
In the rest frame of the quark it is given by T''''(a;)| = 
VAdiag(l,i,i,l)/(127r2ccl) |27]. 
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As discussed in section |lT] the symmetric tensor AT^^ 
is determined by the energy-momentum conservation 
equation (3.11 it obeys together with the perturbations 



in the associated hehcity variables, 



ATo = 



i3 



AT^^^AZ\ , 



?(4) 



(4.58a) 
(4.58b) 
(4.58c) 



where AZf'' , Azf\ and Az''^^ are the vakies of zf^ , 

Z\ and Z\ , respectively, at temperature T minus zero 
temperature. 



E. Numerical Solutions 



One must solve the inhomogeneous differential equa- 



tions for the Z„ Eqs. (4.251, (4.301 and (4.34), with 



appropriate boundary conditions at the horizon and at 
the boundary. To do so, we construct Green's functions 
Gs{u,u') out of homogeneous solutions, 

Gs{u,u') = g<(w<) 5>(u>) / Ws{u') , (4.59) 

with Ws{u) the Wronskian of gf and gf , and convolve 
with the source. 



Zsiu) 



du' Gs{u,u')Ss{u') 



(4.60) 



The appropriate homogeneous solutions to use are dic- 
tated by the boundary conditions. The differential op- 



erators in Eqs. (4.251, (4.301 and (4.341 have singular 
points at u = and u = Uh- The indicial exponents 
at u = are non-negative integers and those a,t u = Uh 
are -iiiuiUh/^- Vanishing of the metric perturbation at 
the boundary requires that gf{Q) = 0, while the require- 
ment that the black hole not radiate ^47^ implies that 
g^{u) ^ (w— Wft)"*'^"''/"' as u ^ Uh- 



1. Tensor Mode 

The indicial exponents of the helicity two differential 
equation (4.25) at it = are and 4. Hence the homo- 



geneous solution 32 niust have the asymptotic behavior 



as u ^ 0. We fix the (arbitrary) overall 



normalization by requiring that lim„^o (") /" ~ 1- 

At zero temperature, the homogeneous solutions to the 
helicity two equation are 



92.T=0m = 2 

q'^—uj 



(4.61a) 
(4.61b) 



where K2 and I2 are modified Bessel functions. The cor- 
responding zero temperature Wronskian is 



M^2,T=o(w) = - 



8u-^ 



(4.62) 



The radial coordinate of zero temperature AdS space 



(using the metric (2.1 1 with Uh ~ 00) runs over the inter- 
val (OjCxd). To subtract zero temperature contributions 
to the stress-energy tensor it will be convenient to map 
this interval onto (0, Uh) with a change of variables 



C(^^)- 



f{u) 



Putting everything together, we have 



Azf^ ^ 



yW2{u) 

9 ^ 



(4.63) 



(4.64) 



+ Jiu) ^^K2[au)V¥^) ^2,T=0(C("))} , 

where J = d^/du ~ {f — uf')/p is the Jacobian for the 
change of variables ( |4.63 1. 



2. Vector Mode 



The indicial exponents of the helicity one differential 



equation (4.30) at u = are and 3. Vanishing of the 



metric perturbation at the boundary thus requires that 
gf{u) ^ V? as u Q. The overall normalization of gf is 
fixed by requiring lim„_^o,9]^(w)/u^ = 1. 

At zero temperature, the homogeneous solutions to the 
helicity one equation are 



2^2 



and their Wronskian is 

W"i,r=o(u) 



2m3 



(4.65a) 
(4.65b) 

(4.66) 



We therefore find 



AZ- 



(3) 



du 



IW^o(w) 

+ j(.)^gH^i., 

2^(m) 



5i(w) (4.67) 
(^(i*)yg^)5i,T=o(C(u))}. 



3. Scalar Mode 

Finally, the exponents of the helicity zero differential 
equation (4.34) at u = are and 4. Hence g^ must 
satisfy g^Xu)^ u** as u — s- 0, and the normalization is 
fixed by requiring liniu^o 5o~(''^)/w^ = 1- 
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At zero temperature, the homogeneous solutions to the 
helicity zero equation are 



8^2 



3o,T=o(^) 



q'^—uj^ 



and their Wronskian is 



0,T-. 



(4.68a) 
(4.68b) 

(4.69) 



Consequently, 



AZ, 



(4) 



lWo(«) 




du ^ So{u) 



+ J{u) 



? 9 



(4.70) 

K2{i{u)^/q^) So,T=o{^{u))] . 
4- Numerics 



FIG. 2: A plot of the window function w{q). The window 
function does not modify the Fourier space data for unq < 20. 



the direction v. The remaining integrals over q and 9 are 
evaluated numerically. 

The integrands (with the factor of q^ from the mea- 
sure) scale like q times an oscillatory function in the large 
q limit. This reflects on the fact that AT^'^{t, x) diverges 



For a given momentum q, the homogeneous solutions 
gf are^ evaluated by numerically integrating^e ap£ro- {-^^ Ty\x-vt\'' in the vicinity of the quark.' To deal with 

this divergence we multiply the Fourier space data by the 
window function 



priate homogeneous differential equations (4.251, (4.30) 



and (4.341 (without sources) outward from the horizon. 
Then the solutions g'f are evaluated by numerically in- 
tegrating the same equations (without sources) inward 
from the boundary. Given these numerically determined 
homogeneous solutions, the perturbations in the helicity 
variables ATg are evaluated by numerically performing 
the radial integrals in Eqs. ( |4.64| , ( |4.67[ ) and ( 4.70| ), and 
inserting the results in Eq. (4.581. 



w{q) 



rf 



/ Uhq - 30 \ 
V 4.5 / 



(4.74) 



The reconstruction of the stress-energy tensor compo- 
nents from the helicity variables is given in Eqs. (2.8) 



and (2.9). In particular, the perturbation in the energy 



density is 



A£{t,x) 



(2^ 



ATo{uj,q) e 



iQ-x 



(4.71) 



while the perturbation in the energy flux is 



AS{t,x) 



d^q 
(2^ 

iQ-x 



;ATo-z^°)+Afi(c^,q) 



X e 



(4.72) 



The frequency integrals in these expressions are trivial 
as the integrands are proportional to 2tt5{uj—v ■ q). We 
write the remaining 3(i integrals in spherical coordinates 
(g, 6, (p) with 9 the polar angle measured relative to the 
direction of the velocity v. In this coordinate system we 
have 



qcost 



q± = q sm f 



(4.73) 



This window function is shown in Fig.[2j From the figure 
we see that the window function does not significantly 
modify the Fourier space data for Uhq < 20. Corre- 
spondingly the window function alters the real space en- 
ergy density and energy fiux over length scales Uh/20, 
which is much smaller than the position space grids used 
to make the plots shown below. With this window func- 
tion, the integral over q is evaluated with an upper limit 
of Uhq — 40. 



V. RESULTS 

For small distances d = \x—vt\ <^1/T away from the 
moving quark, the dominant contributions to the stress- 
energy tensor come from the T = stress-energy tensor, 
which scales like 1/c?^. To highlight the medium depen- 
dent perturbations in the stress-energy tensor we defined 
AT^^'^ [in Eq. (4.56)] to be the temperature dependent 
perturbation to T^"^. 

Figures [3j|5] show plots of AE{t, x) and AS{t, x) = 
\AS{t,x)\ at quark velocities v = 1/4, v = 3/4 and v = 
l/\/3, respectively.^^ The flow lines superimposed on the 



where q\\ is the component of q parallel to the quark's ve- 
locity. The integration over (j) is easily done analytically, 
as the physical problem is cylindrically symmetric about 



In making these plots, we use spatial grids with resolution Ax 
equal to one to two times 1/(2itT). This limits the fidelity of 
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|x|Ag(x) ! i \ |x|A5(x) 




FIG. 3: Left— Position space plot of \x\A£{x)/{T^^/X) for v = 1/4. Right— Position space plot of \x\AS{x)/{T^^/\) for 
V — 1/4. The flow lines on the surface are the flow lines of the energy flux AS{x). There is a surplus of energy in front of the 
quark and a deficit behind it. Correspondingly, trailing the quark there is a stream of energy flux which moves in the same 
direction as the quark. Note the absence of structure in A£{x) for distances \x\ 3> l/(7rT). 



|x|A£(x) 




FIG. 4: Left— Plot of \x\A£{x)/{T^^/X) for v = 3/4. Right— Plot of \x\AS{x)/{T^y^) for v = 3/4. The flow lines on the 
surface are the flow lines of AS{x). There is a surplus of energy in front of the quark and a deficit behind it. Correspondingly, 
trailing the quark there is a narrow stream of energy flux which moves in the same direction as the quark. A Mach cone, with 
opening half angle 9m ~ 50° is clearly visible in both the energy density and the energy flux. Near the Mach cone, the bulk of 
the energy flux flow is orthogonal to the wavefront. 
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FIG. 5: Left— Plot of \x\A£{x)/ [T^y/X) for v = c^. Right— Plot of \x\/\S{x) / {T^ y/X) for v = c^. The flow lines on the surface 
are the flow lines of the energy flux AS{x). A planar Mach cone is visible in both the energy density and the energy flux. Near 
the Mach cone, the bulk of the energy flux flow is orthogonal to the wavefront. 



plots of S{t^ x) are the flow lines of S{t, x). In these plots 
the quark, at the time shown, is at a; = 0. We use units in 
which ttT = 1. Since A/" = 4 SYM is a conformal theory, 
the speed of sound is l/-\/3- Hence Fig.jsjshows subsonic 
motion. Fig. [4] shows supersonic motion, and Fig. |5] is 
precisely at the speed of sound. 

As discussed earlier, AT^'^ may be reconstructed from 
the helicity variables A7^, combined with the energy- 
momentum conservation equation (3.1) and the vanish- 
ing trace condition of the SYM stress-energy tensor. In 
Appendix IS] we compute the small momentum limit of 
ATs to 0(q^ with the ratio r = uj/q held fixed. Defining 



{2^)5{uj-vq), (5.1) 



one finds 
ATo = 3/C 



27rii^\/l — 



(l-3r2)q (l-3r2)2 
vqj_ \ J_ _^ Uft(l-4r2)" 
irq 



4^2 



AT, = -/C ^ (ei ® ei - ® e,), 

2 g-^ 



(5.2a) 
(5.2b) 
(5.2c) 



up to 0{q) corrections. 

These results for the long wavelength limit of ATg 
may be compared to the analogous long wavelength 



these plots for distances \x\ ~ Ax from the quark. As noted 
earlier, the temperature dependent energy density and energy 
flux behave like T^/\Ax\'^ for distances |Aa;| < 1/T from the 
quark. See Refs. |27l 149) for a discussion of the stress-energy 
tensor in the near zone. 

These results for the small momentum limit of the stress-energy 
tensor may also be found in Ref. |27| . 



limit of the hydrodynamic quantities AT^''^'^''" given in 
Eqs. (3.231. To give a meaningful comparison, we ex- 



pand the hydrodynamic helicity variables in powers of q 
with the ratio r = uj/q fixed. Writing 



P - 
Jt 



qp 



(1) 



(0) 



-qji^U 



(5.3a) 
(5.3b) 



we find 



. ^hydro ^ 3p(i) _ 9»r7pW + 3(l-3r2)p(2) 
° (l-3r2)g (l-3r2)2 



r(0) n /-(") I ,> r(i) 

hydro "^T , ^«'T ' "^T 



ATI 



irq 



(5.4a) 
(5.4b) 
(5.4c) 



neglecting 0{q) corrections in Eqs. (5.4a) and ( 5.4b[ ) 



hydro 



to 



Comparing the expansions of ATq and A7^ 
that of ATo and AT, we see they agree provided the 
transverse momentum diffusion constant and the sound 
attenuation constant have the expected values given ear- 
lier,2" namely D = l/(47rr) and 7 = l/(37rr), and the 
hydrodynamic sources have the expansions 

p = /C [ iuj{l + u^) _ UhUj'^ + 0{q^)\ , (5.5a) 

vq±_ 



Jt = /C ^ [1 iuuLo + 0{q^)\ ii 



(5.5b) 



Higher order transport coeflicients such as 0, the coefficient mul- 
tiplying second derivatives of the energy density, can be deter- 
mined by carrying out the expansion of the helicity variables to 
higher order in q. 
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1 dPn 




30 60 90 120 150 180 

e 



FIG. 6: Plot of {T'^^/xy^dPR/d.e for v = 1/4 at distances R = 5/(7rr) (black curve and red triangles), 10/(7rT) (blue curve 
and purple squares), and 2Q/{ttT) (green curve and black dots). Solid lines are the full AdS/CFT results, while the symbols 
show the linearized hydrodynamic approximation. Note the presence of energy being radiated into the forward hemisphere 
in front of the quark. Negative values in the backward hemisphere are due to the diffusion wake. At all three distances, the 
hydrodynamic approximation to dPn/dO agrees very well with the AdS/CFT results. 



Using the fact that J° = F", Eqs. (3.241 and (5.5a) imply 



that the longitudinal component of the effective source 
has the expansion 



q 



Jl = /C ^ [tJ + iuhOj"^ - ijv^q^ + 0{q^)] 



Combining Eqs. (5.5b) and (5.6), we find 



J = IC [v + iuhLov 



(5.6) 



(5.7) 



Comparing with Eqs. ( 3.2 1 and ( 3.3 1 , one sees that, as ex- 



pected, the leading term in the gradient expansion of the 
effective source J is simply the microscopic force density 
F. When Fourier transformed to position space, sublead- 
ing corrections to the effective source generate derivatives 
of delta functions at the location of the quark. 

As the above makes clear, the small q expansions of 
ATq and ATi are completely consistent with their hy- 
drodynamic counterparts. However the small q limit of 
AT'^'^, as reconstructed with the helicity variables ATg, 
does not agree with the small q expansion of ATj^^^^.^. 
There are two sources of discrepancy. First, the helic- 
ity two variable AT2 is non-vanishing at 0{q'^) while 
^^hydro yg^-^jg]-^gg identically in the large limit. Sec- 
ond, even with identical helicity variables, the reconstruc- 
tion of AT^'' from its helicity variables ATs differs from 
the reconstruction of AT^^^^^ from AT^'^^'*''", because 
AT^"^ and ATl^_^^^.^ satisfy differing energy-momentum 
conservation relations. In particular, AT^'' satisfies the 
exact microscopic conservation relation (3.1) with the 



force density F" on the right hand side, while AT^ 



satisfies the effective energy-momentum conservation re- 
lation (3.12) involving the effective source J'^. However, 
examining the small q limits of AT^'^ and AT^^^^^, we 
find 



where 



^ hydro ^ ' 

A'^ = IC (v'^-fSij - Uh ViVj) 



(5.8) 

(5.9a) 
(5.9b) 



When Fourier transformed back to position space, A^'^ 
has point support at the location of the quark and satis- 
fies 



(5.10) 



In other words, A'^'^ is a local contribution to the stress- 
energy tensor which is precisely tailored to compensate 
for the difference between the effective source and the 
microscopic force density. This is a necessary result. The 
gradient expansions of AT^^" and AT{^^^^^ must agree far 
from the quark. However, in the long wavelength limit, 
any discrepancy between AT'^'^ and AT-^^^^^ in the near 
zone may be represented as a term with point support 
at the location of the quark. Because AT'^" and ATj^/^^.^ 
satisfy different conservation equations, A'^'^ must obey 
Eq. (5.10) to all orders in the gradient expansion of J^. 



hydro 



It is, of course, reassuring to see A'^'^ naturally emerge 
from the gravitational calculation. Note that informa- 
tion from all three gauge invariants Zg is used to calcu- 
late A'^'^ , thus showing the non-trivial interplay between 
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FIG. 7: Plot of {T^y/\)-^dPR/de aX v = Cs and distances R = 5/(7rr), 10/(7rr), and 20/{ttT). The labeling of curves is the 
same as m Fig. [6] Note the presence of energy being radiated into the forward hemisphere in front of the quark, together with 
the inward flux associated with the diffusion wake in the backward hemisphere. The Mach cone is centered at ^ = 90°. At all 
three distances, the hydrodynamic approximation to dPR/d9 agrees very well with the AdS/CFT results. 



1 dPg 

T^Vx de 



-3- 



— R^5/(ttT) 
A Hydro at R = 5/(ttT) 

— fl = 10/(7rT) 
= Hydro at i? = 10/(7rr) 

— R^20/(7tT) 
. Hydro at R = 2Q/{itT) 
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FIG. 8: Plot of {T^VX)-^dPR/de at v ^ 3/4 and distances 7? = 5/{nT), 10/(7rT), and 20/{nT). The labeling of curves is the 
same as in Fig. [6] Note the absence of energy radiated in front of the quark, as the quark is moving supersonically. The Mach 
cone near 130°is sharper than in the transonic case shown in Fig.[7] The hydrodynamic approximation to dPn/dO improves as 
R grows, but is already fairly good at i? = 5/{nT). 
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the decoupled gravitational equations of motion (4.251 
p30|, and (OSl). 



It is illuminating to compare hydrodynamics to the 
AdS/CFT results in position space. A simple quantity 
to compare, which involves both the sound and diffusion 
modes, is the instantaneous angular distribution of the 
energy flux at a distance R from the quark. Choosing to 
evaluate this at t — 0, when the quark is at the origin, 
the angular distribution of the energy flux is 



dPR 

dn 



R'^x-AS{0,RS:). 



(5.11) 



Because of the rotational symmetry about the velocity 
axis, its sufficient to consider the energy flux per unit 
polar angle 9, 



dPR 
dO 



= 2ttR^ sin 9 (cos 9 AS« + sin 9 AS±) 



(5.12) 



with 9 = corresponding to the v axis and AS\\ and 
AS'j^ the components of the energy flux parallel and per- 
pendicular to the velocity vector v, respectively. We plot 
dPii/d9 for V = 1/4 in Fig. [6] w = Cg in Fig. [t] and 
V = 3/4 in Fig. |8] In all three flgures we plot the an- 
gular distribution of the flux at distances i? = 5/ (ttT) , 
10/(7rr), and 20/(7rT). The hydrodynamics curves were 
made by numerically integrating Eqs. (3.271 and (3.28) 



with the leading order effective source J'' = F''. Hydro- 
dynamics becomes increasingly accurate at longer dis- 
tances, but even at i? = 5/{nT) one sees rather good 
agreement between linearized hydrodynamics and the full 
AdS/CFT results. 

We emphasize that dPj^/dil is not the same as the 
power radiated in the rest frame of the quark. In the 
quark's rest frame, the total power radiated through 
spheres of radius R is independent of R. However, in the 
rest frame of the plasma, the total flux, J dil {dPa/dn), 
through a sphere centered on the instantaneous posi- 
tion of the quark grows with increasing R as successively 
larger spheres capture energy radiated ever farther back 
in time. 



VI. DISCUSSION 

Much of the qualitative and quantitative structure in 
the plots of A£{t,x) and AS{t,x) can be understood 
from hydrodynamic considerations alone. This statement 
is reinforced by Eq. 



(5. 



which shows that the long 
wavelength limit of the stress tensor as computed with 
gauge/string duality coincides with the linearized hydro- 
dynamics result, provided the latter is computed with 
the correct effective source. As discussed in Section [HT] 
long wavelength perturbations in the energy density sat- 
isfy the diffusive wave equation (3.291 describing sound 
waves with the dispersion relation 
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FIG. 9: A schematic representation of the Green's function 
solution to the sound equations (3.291 and (3.351 for super- 
sonic motion. As the quark moves, it creates sound distur- 
bances which propagate out in spherical shells. The result- 
ing sound waves add coherently along the Mach cone, shown 
in blue. Because the sound waves are damped, each spher- 
ical wave will broaden as it propagates. The broadening is 
schematically indicated by the dashed lines. This implies that 
the Mach cone broadens with increasing distance from the 
quark, as indicated by the red lines. 



[with Cs = l/\/3 and 7 = l/(37rr)]. A textbook con- 
structive interference argument, illustrated in Fig. [9] 
shows that a projectile moving supersonically will pro- 
duce a Mach cone with an opening half-angle given by 
sin^M = Csjv (where tan 0m = For u 3/4 

this is 9yy = 50.3°. As is evident in Fig. [i] at w = 3/4 
the energy wake is concentrated, as expected, along a 50° 
cone with the the bulk of the associated energy flux flow- 
ing perpendicular to the wave front. Similarly, in Fig. [5] 
at w = Cg the energy wake is concentrated along the plane 
X|| = with the bulk of the associated energy flux again 
flowing perpendicular to the wave front. 

The damping of sound waves implies that their wave- 
forms must broaden as they propagate. This in turn 
implies that the Mach cone broadens with increasing dis- 
tance from the quark. This behavior is clearly seen in 
Fig. |4] As one can reason from Fig. |9] the width FMach 
of the Mach cone, defined as the length scale over which 
the energy density attenuates (exponentially) in the for- 
ward direction, increases with distance d from the quark 
like 



Mach 



(7d) 



1/2 



(„2_c2)l/4 



(6.2) 



±Csq 



79 



(6.1) 



This attenuation length diverges as v Cg', in this limit 
there is only power law falloff of the Mach cone in front 
of the quark. This reflects the fact that at v = Cs, sound 
waves emitted by the quark arbitrarily far in the past 
can add coherently to the Mach cone at finite distances 
from the quark. If a given sound wave has to travel a 
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this more formally note, from Eq. (3.34b I, that the diver 



FIG. 10: A schematic representation of the Greens func- 
tion solution to the diffusion equation ( |3.37[ ). As the quark 
moves, momentum deposited in the diffusion channel at time 
to diffuses over the length scale \/D(t — to). This results in 
the broadening of the diffusion wake with increasing distance 
from the quark. 



distance £ to reach the Mach cone, the width of the wave 
will have broadened to ~ y/ji/cs- At precisely v — Cg, 
sound waves with i oo contribute to the Mach cone 
at finite distances from the quark, leading to power-law 
falloff of the Mach cone. 



From Eq. (3.371, we see that part of the long wave- 



length perturbation in the energy flux satisfies a dif- 
fusion equation with diffusion constant D = l/(47rT). 
Fig. [To] shows a pictorial representation of the the G reens 
function solution to the diffusion equation (3.371. As 



the quark moves, momentum transferred to the diffusion 
mode at time t mode will gradually diffuse outward, re- 
sulting in a stream of momentum trailing the quark. The 
width of the corresponding diffusion wake grows with dis- 
tance d from the quark like 



ffusioii 



y/Dd/v. 



(6.3) 



This behavior is evident in the energy flux plots in 
Figs. [3}|5j The directional dependence of the diffusion 
wake is dictated (in the long wavelength limit) by the 
directional dependence of the effective source Jdiffusion, 
(defined in Eq. 3.381. To leading order in gradients, 
•^diffusion oc v, SO the diffusion wake should flow in the 
same direction as the quark's motion. Again, this behav- 
ior is clearly seen in the energy flux plots of Figs. |3}]5] 

It is noteworthy that the diffusion wake is not seen at 
all in the plots of the energy density. This can be un- 
derstood from hydrodynamic and large Nc scaling con- 
siderations. First of all, note that the kinetic energy 
associated with bulk motion of the fluid is negligible 
in the large Nc limit. This follows from the fact that 
the fluid velocity which results from the motion of the 
quark through the plasma is 0{1/N^). (The fluid ve- 
locity u equals the O^N^) momentum density divided 
by the 0{Nc) enthalpy, £+p-) And hence the associated 
bulk kinetic energy density, given by (e+p) , is also 
0{1/N^).'^^ This means that an 0{N°) momentum den- 
sity need not produce any imprint on the O^N^) energy 
density perturbation, exactly as seen in Figs.[3]45] To see 



Consequently, in the large Nc limit, all of the energy lost by a 
quark traversing the plasma is turned into heat. 



gence of the diffusive energy flux has point support at the 
quark. That is, V • 5'diffusion — everywhere away from 
the quark. This implies that the momentum transported 
via the diffusion flux does not involve any corresponding 
perturbation in the energy density. This is analogous, 
in electromagnetism, to the possibility of having a trans- 
verse electric current density without any corresponding 
charge density. 

To put our discussion on a somewhat more quantita- 
tive footing, we now turn to the plots of dPji/dO shown 
in Figs. [6||8] For subsonic and transonic motion (Figs. [6] 
and|7]), dPn/dO is reproduced quite well by hydrodynam- 
ics all the way down to distances R = 5/{ttT) from the 
quark. For supersonic motion (Fig. [8| the agreement 
between hydrodynamics and the complete result is not 
quite as good, quantitatively, at i? = 5/{ttT), but the 
level of agreement is still rather remarkable. The region 
of largest discrepancy occurs on or near the Mach cone. 
It is natural that the hydrodynamic approximation to the 
stress tensor deviates more for supersonic motion than for 
subsonic motion. This is because gradients in the energy 
density and energy flux are largest on the Mach cone. 
This behavior is clearly seen in the near zone behavior 
of the energy density in Fig. [4] However, in constructing 
the hydrodynamic approximation to the stress tensor in 



Eq. (3.11), we have only kept the leading viscous terms 



in the gradient expansion and neglected Q and all other 
higher order terms in the gradient expansion of the stress. 

We emphasize the importance of including the non- 
zero viscosity when evaluating the hydrodynamic curves 
shown in Figs.[6}|8j Had we neglected viscosity, the stress 
tensor would be discontinuous on the Mach cone and, 
for V — Cg, would actually diverge on the Mach cone. 
Furthermore, in the limit of zero viscosity the width of 
the diffusion wake vanishes, leaving a diffusion wake pro- 
portional to S'^{x±)Q{—xi^\ + vt). Therefore, neglecting 
viscosity yields a very poor approximation to the stress 
tensor in the vicinity of both the Mach cone and the dif- 
fusion wake. 

The agreement between hydrodynamics and AdS/CFT 
at distances d>l/T from the quark should be contrasted 
with the corresponding situation at weak coupling. If the 
't Hooft coupling A <C 1, then the SYM plasma has a 
quasi-particle description. The mean free path of quasi- 
particles (gluons, fermions, or scalars) scales like [33] 



1 



f-infp 



TA^lnA-i 



(6.4) 



Therefore, the mean free path, at weak coupling, is para- 
metrically longer than 1 /T. Hydrodynamics is only valid 
on spatial scales large compared to the mean free path. 
So in contrast to the strong coupling results discussed 
above, in a weakly coupled plasma hydrodynamics is 
never valid on distance scales comparable to 1/T. 

It is interesting to ask how the energy and momentum 
lost by the quark is distributed into the sound and dif- 
fusion modes. This question has been addressed in the 
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steady state limit in Refs |31l I32j . However, for reasons 
which are outlined below, we question the physical value 
of addressing this topic in the limit where the quark has 
been moving forever. 

Let V denote a large sphere of radius R whose location 
is fixed and centered on the position of the quark at time 
t = 0. The rate of change of the energy inside volume V 
is 



dEy 

dt 



(6.5) 



After using the energy-momentum conservation equation 
(3.1 ), and separating the energy flux into sound and dif- 



fusive pieces, Eq. (6.5) becomes 



dEv 
dt 



dY.-S, 



sound 



dV 



ffusion 



dV 



(6.6) 



The first term on the right is simply the rate at which the 
quark deposits energy into the plasma while the second 
and third terms are the rates at which energy is added to 
and removed from V via the sound and diffusion fluxes. 

In the steady state limit, it is easy to calculate the rate 
at which energy is added to V by the diffusion flux. If 
R 3> 1/r, then one may compute the surface integrals 
in Eq. (6.6) using hydrodynamics. The long wavelength 



limit of the diffusion flux satisfies^^ 



diffusion 



(6.7) 



From this expression we see that the diffusion flux adds 
energy to the volume V at the rate 



• S'diffusion — 



av 



1,2 • 



(6.8) 



The source of the energy supplied to V via the diffusion 
flux comes from the quark itself — in the distant past 
when the quark was outside of the volume V , it created 
the diffusion flux which subsequently flows into V . (The 
fact that the diffusion flux adds energy to V simply re- 
flects the fact that the diffusion flux flows in the same 
direction as the quark's velocity.) Note that the rate at 
which energy is added to V via the diffusion flux is inde- 
pendent of the size of V , and so is non-vanishing in the 
R oo limit. Furthermore the rate (6.i 
[recall that /" = 



oo limit 
non-zero limit as v 



has a flnite. 



So the diffusion flux adds energy to V even when the 
quark's velocity is taken to be arbitrarily small! The ori- 
gin of this peculiar behavior comes from the assumption 



that the quark has been moving forever. In this limit, 
diffusive energy flux deposited in the arbitrarily distant 
past can influence the rate of change of the energy in the 
volume V. Completely analogous conclusions also hold 
for the sound flux.^'^ 

The assumption that the quark has been moving at 
constant velocity forever presents a great technical sim- 
pliflcation in both the gravitational and hydrodynamic 
computation of the quark wake. However, the utility of 
this simplifying assumption is limited to questions which 
are insensitive to the details of the quark's trajectory in 
the distant past. Asking how the energy and momen- 
tum lost by the quark are distributed into the sound and 
diffusion modes, when the quark has been moving for- 
ever, is not such a question. A better and still analyti- 
cally tractable question is to consider is how much energy 
and momentum are deposited into each mode when the 
quark has been moving at constant velocity v for a long 
but flnite period of time At — tf — ti. The total energy 
transfered to the plasma is then 



AE{t) = j d^x 



{t,x) 



(6.9) 



where the integration is taken over all of space. For late 
times after the quark's motion has ceased, all of the en- 
ergy deposited by the quark will be transported away via 
sound waves. That is, because the quark has been mov- 
ing for a flnite period of time, no fluxes at spatial inflnity 
contribute to Ai?(i). 

The total momentum deposited by the quark is 



Ap = Atf. 



(6.10) 



To see what fraction of the quark's momentum is de- 
posited in each mode, we will examine the diffusion mode 
and then infer the momentum deposited in the sound 
mode via momentum conservation. The total momen- 
tum transfered to the diffusion mode is 



Apdiffusion(t) = J d^X S'diffusion(t, x) 



(6.11) 



If the interval during which the quark has been moving 
is sufficiently long, then the total momentum transfered 
to the diffusion mode may be computed with hydrody- 
namics. Moreover, in the limit At ^ l/T, the vast ma- 
jority of the volume of the diffusion wake will be well ap- 
proximated by its steady state limit. In the steady state 
limit it is easy to solve the diffusion equation (3.37) with 



This equation was first obtained in Refs. 1311 I32| via a grav- 
itational calculation. However, it follows directly from lin- 
earized hydrodynamics using the leading order effective source 
JM = piJ- , Specifically, Eq. (|6.7[l follows from the definition 
of the diffusion flux in Eq. | |3.34b[ l plus the observation that 
V • Snoniocal = — / , to leading order in gradients. 



To demonstrate this, one may solve the wave equation | |3.35| in 
the analytically tractable limit of small velocity. When v <^ Cs 
no Mach cone is generated and the effects of viscosity on the 
long distance behavior of the sound energy flux are negligible. 
Taking the duration of time in which the quark has been moving 
At — > oo before taking ?; — > 0, one finds that the rate at which 
the sound flux removes energy from the volume V has a finite, 
non-zero limiting value. 
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the corresponding leading order effective source given in 
Eq. ( |3.39b I . Assuming that the quark is at a; = at time 
tf, the result reads 



'S'difFusion(i — tf , x) 



AnDr 



exp 



2D 



(^11 + ''") 



(6.12) 



where r = \x\ and the ellipsis denotes corrections sup- 
pressed by inverse powers of r (coming from higher order 
gradient corrections). Eq. (6.12) provides a good approx- 



imation to the diffusion wake everywhere except in a re- 
gion of size ^ V DAt centered about the quark's starting 
point. In the At — > oo limit this region is negligible com- 
pared to the length of the diffusion wake, which is v At. 
To obtain the total momentum transfered to the diffu- 
sion wake, we integrate Eq. ( 6.12| over a length v At in 
the a;|| direction. Specifically, we have 



course, the application of hydrodynamics requires among 
other things, knowledge of the viscosity of the plasma and 
knowledge of the drag force on the quark, neither of which 
are under very good control for real quark-gluon plasma 
at accessible temperatures. (And, as noted earlier, ne- 
glecting viscosity will yield a poor approximation in the 
vicinity of both the Mach cone and the diffusion wake.) 
However, at least when Nc is large (and the hydrody- 
namic equations are linear) , the structure of quark's wake 
is rather insensitive to magnitude of the drag force. For 
example, in the steady state limit the magnitude of the 
drag force enters only as an overall normalization of the 
quark wake. As a next step toward a more realistic treat- 
ment, it would be interesting, and should be feasible, to 
study the wake produced by a quark traversing an ex- 
panding and cooling A/" = 4 SYM plasma [16j . 



Apdiffusion(t/) = / dx\\ d X± Sdmusion{tf,x) . 

J-vAt J 

(6.13) 

The integral is elementary to carry out in the At — > oo 
limit. The result reads 



Ap. 



'diffusion 



itf) = Atf. 



(6.14) 



This shows that in the At — > oo limit all of the mo- 
mentum lost by the quark is transfered to the diffusion 
mode. Correspondingly, the sound mode carries no net 
momentum. The qualitative origin of this statement can 
be understood by inspecting the momentum density near 
the Mach cone in Figs. |4]-|5] In particular, the momen- 
tum density (i.e., energy flux) flows predominately paral- 
lel and anti-parallel to the Mach cone normal. Evidently, 
when integrated over all space, the opposing flows cancel. 

The direct applicability to RHIC physics of results ob- 
tained from studying the wake produced by a quark mov- 
ing for an unboundedly long time through a plasma of 
infinite extent is, of course, questionable at best. The 
plasma produced at RHIC is very dynamic. The expan- 
sion of the plasma and the resulting variable speed of 
sound can have a dramatic effect on the structure of the 
quark wake |24j . However, the agreement between hy- 
drodynamics and gauge/string duality in describing the 
structure of the quark wake in strongly coupled SYM 
does have relevance for RHIC. In particular, at least for 
the particular problem addressed in this paper, it shows 
that one can use hydrodynamics to address physics at 
length scales all the way down to distances less than two 
times 1/T. 

Because the quark-gluon plasma produced in heavy 
ion collisions is beheved to be strongly coupled [TJ [5], 
the agreement between hydrodynamics and gauge/string 
duality in a strongly coupled SYM plasma bolsters the 
notion that one should be able to model accurately the 
wake produced by a high energy quark (or gluon) moving 
through a QCD plasma merely using hydrodynamics. Of 



VII. CONCLUSIONS 

Using gauge/string duality, we have evaluated the per- 
turbation in the stress-energy tensor due to the pres- 
ence of a fundamental quark moving through a strongly 
coupled large maximally supersymmetric Yang-Mills 
plasma. Our plots of the energy density and energy flux 
in Figs. |3}]5] clearly display the formation of Mach cone 
for velocities v > Cs, together with the presence of a dif- 
fusive wake in the energy flux at all velocities. 

We have argued that the effective source for hydro- 
dynamics is, to leading order in gradients, simply deter- 
mined by the microscopic drag force acting on the quark. 
By comparing the small momentum asymptotics of the 
stress-energy tensor with the predictions of hydrodynam- 
ics, we were able to determine the flrst subleading cor- 
rection to the effective hydrodynamic source. 

Using the leading order effective source for hydrody- 
namics, we compared the hydrodynamic prediction to 
the complete result obtained via gauge/string duality for 
the instantaneous angular distribution of power radiated 
through spherical shells a distance R from the quark. The 
comparison showed remarkably good agreement between 
hydrodynamics and the exact result at distances all the 
way down to i? ~ 1.6/r from the quark. This reinforces 
the notion that wakes produced by high energy particles 
traversing a real quark-gluon plasma can be well modeled 
merely using hydrodynamics. 
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APPENDIX A: THE BOUNDARY ACTION 

The total gravitational action for the heavy quark ef- 
fective theory is given by 



5, 



EH 



Sr 



GH + '5'dBI + -^NG + "^CT 



Sk 



Sc 



where 



5eh 



1 
1 
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d^x V-G {R+2A) 



d^x V-7 2K , 



(Al) 

(A2a) 
(A2b) 



and the other pieces were described in section IV We 
write the metric in the bulk as 



Gain — G^^j^ + Hmn , 



(A3) 



and choose to work in a gauge where h^M = 0. We 
similarly write the metric on the boundary as 



hanging down to the horizon. The variation of the 
Nanibu-Goto action is, by definition 



-0(0) ^Sh 



MN 



MN 



(All) 



Hence the variation in the Nanibu-Goto action cancels 



the last term in Eq. (A9l. As discussed in section IV 



in the large quark mass limit, the DBI action may be 
replaced with the ordinary Maxwell action for the elec- 
tromagnetic field (residing on the boundary) which ac- 
celerates the quark. 

The trace of the extrinsic curvature is given by 



K = Wnn 



M 



(A12) 



where n*^ — —\/G^ 5^^'^ is the outward pointing normal 
to the boundary and denotes covariant differentiation 
with respect to the full metric G^^. Using 



^7 Vun 



M 



(A13) 



To quadratic order in we have 



-G(i? + 2A) = ./-G(o) iDpW'^ + Ca- 



2d\ 



(A4) 



(A5) 



where 



Co ^ iDMhO^'h+lDphMND^^h 



MiNP 



iDph^'^DMh - \Dph'''''DPhMN 



2L2 



uMN, \ 



(A6) 



and 



\hDPh+\DM{hh^''') 

},hP''DMh + h''''DPhMN 



n-MN) 



- D^' {h 

To vary the gravitational action, let 



(A8) 



where /i^jy satisfies the linearized equations of motion 
(4.181 and bh^j^y is infinitesimal. The variation of the 
bulk action density is 



— G(o) (5/3o 



where 



= \8hD''h - \bhMN D^h^'^ + ShMN D^'hP'' 

-\6h DnH^^ - \5hP^ DnH , (AlO) 

and t*^^ is the bulk stress-energy tensor, whose only 
source (in the limit of large quark mass) is the string 



we can write 



S-GH = -4 lim \/G55 du I d'^x 

4 M^O J 



-7- 



(A14) 



Assembling the pieces, and specializing to the AdS- 
Schwarzschild bulk geometry, we have 



5Sq = 5Sb, + 5S] 



horizon ; (-^-^^) 



where 



5Sv 



lim 

24 "^0 



d^x 



III 
2u3 



7-3 f 



xSh^,,h^''-^6hh-^dhooh- 
2 Afu 



2u3 



Shot K 



3 (V/ -l)L 



dhn 



xSh' 



(55, 



EM : 



(A7) with 



(-6/ + 677 

21*4 



(A16) 



(A17) 



and ^S'horizon is a surface term at the black-brane hori- 
zon whose explicit form will not be needed. From the 
linearized field equations (4.181 it is easy to see that near 
x^. Using Sh^i, = 



the boundary h 

that all of the terms in the curly braces in Eq. (|A16| scale 



SH^i^, we see 



like u SHfj^i, as u ^ 0. It follows that these terms will 



(A9) ^'^^ contribute to SSq/SH^^^x^u) in the u 



Neglecting terms which vanish in the u 
therefore have 



limit. 
limit we 



(55b — lini 



d'^x 



4K5U3 



(55, 



EM ■ 



(A18) 
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where jy^'' = diag(— 1, 1,1,1) is the Minkowski space met- 
ric and 

T,';'^ = ediag(l,i,|,|) (A19) 

is the stress-energy tensor of the equihbrium A/" = 4 SYM 
plasma. The expHcit value of the equilibrium energy den- 
sity is 



3L^ 
2^ 



(A20) 



We write the variation of the boundary action in 
Fourier space and wish to express it in terms of the gauge 
invariants Zs- The variation in the boundary action is 



(5S'r = lim 



n^O ] (27r)4 



J__ ( _ 



EM 



(A21) 



Here products are understood to be evaluated at an- 
tipodal momentum as dictated by the reality condition 



Tensor Mode 



The terms in Eq. (A21I constructed out of radial 
derivatives of Zf'' = Hat — 2^ab'H are 



where 



A. 



A2 SZfduZf, 



(A22) 



(A23) 



Vector Mode 



3. Scalar Mode 



The terms in Eq. (A21I constructed out of radial 
derivatives of Hoa , Ti-aa , Ti-qq and Hoq are 



5S?> — lim f . ^, „ „ „ 
^ u^o J (27r)4 8nlu^ 



25Hqq duTiqq — STiaa dvHaa 



(A28) 



Using the first order equations (4.33b |-( 4.33d I, this is 
seen to be 



^S'r = lim / - — ^ 
^ u^o J (27r)4 



00 



where 



and 



^SHqq d'" + \5naa T'^ + in^q 

L3 



lOO 



Aq = 



-3lu 



(5q2-3w2) to 



q {q^-'Suj"^) tq5 + iuq^ (q'^-uj^) t^^, 



UJ iuj'^-'iq^) % 



q (5w^-3g^) tg5 + iuuj^{q'^-uj^) t^^ 



iu{q^ —uJ^) t 



^05 + 9 tqb 
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q iuj -3q ) to 



■ LO (g^-3w^) ig5 -I- iuq {q^~u?) t^^ 



(A29) 
(A30) 

(A3 la) 
(A31b) 
(A31c) 
(A31d) 



The terms in Eq. (A21l constructed out of radial 



derivatives of Tiaq and Tioa are 

(A24) 

We now use Eq. (4.291 to eliminate duHqa- Doing so, we 
find 

6Si = liin I 7;^[A (qSHoa + C^SHaq) Zl+SUaqd'"'] , 

(A25) 
(A26) 



u^o J (27r)4 



where 



Ax = 



and 



7-3 

= ^ ta5 . 



{A27) 



APPENDIX B: SMALL MOMENTUM 
ASYMPTOTICS 

For very large or very small momentum (compared to 
T), one may find explicit asymptotic expressions for the 
homogeneous solutions gf and gf , and derive the result- 
ing asymptotic behavior of the SYM stress-energy tensor. 
In what follows we only consider the small momentum 
limit and, in particular, calculate Tg to C(l). 

We write the homogeneous solutions gf and 5^ as a 
power series in q, 

-fg0W + g2^(2)+O(g3), 



(Bla) 
(Bib) 



The differential equations (4.251, (4.30) and (4.341, with- 



out sources, are then solved order by order in q subject 
to the boundary conditions discussed in Section |IV E| In 
what follows we define r = uj/q. 
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1. Tensor Mode 



and 



For the tensor mode it is s ufficie nt to expand the ho- 
mogeneous solutions to Eq. (4.251 to 0{q'^). Doing so 
one obtains 



4°) = -4 In/, 



and 



(0) 



1. 



(B2) 



(B3) 



The corresponding small momentum limit of the Wron- 
skian is 



W2 = -—+0{q). 



(B4) 



Substituting the ab ove expansions into the definition 
of AZj^"* in Eq. (4.641 and performing the radial integral, 
one finds 



AZ. 



(4) 



ivq_ 



X {ii ii - 62 £2) + 0{q) . 



(B5) 



Substituting AZ^^^ into Eq. (|4.54c| yields 



AT, 



^ {2n)d{uj - V ■ q) 



X (ei (g) ei - €2 «i £2) + 0{q) . 
2. Vector Mode 



(B6) 



(1) 



iru iru 



2< 



(2+^ln^) 

2u^ , _i u2 
— tanh 



(B9a) 
(B9b) 



u\ f 



u 



44- In^ -8tanh"^ 



u 



16 f f 

In21n4-ln— f-ln/_ +lnl61n^ 



+ 21n(//_)ln/+-lnV++4Li2 
The corresponding Wronskian is 

,3 



. f- 



Wi 



7 



irq 



Uh 



(l + r2ln4) q^ + 0{q^) 



(B9c) 



(BIO) 



Substituting the above solutions into Eq (|4.67|, one finds 

1 



AZ. 



7(3) 



4V\ 



vqj_ 



irquh 



4^2 



X {2tt)5{u} — V ■ q) ii . 
From Eq. ( |4.54b I we therefore have 



ATi 



V A vqj^ 
2TTUhVl-v^ q irquh 
X (27r)(5(cj ~ V ■ q) ei 



l_4r2 



+ Oiq) 
(Bll) 

Oiq) 
(B12) 



For the vector mode it is necess ary to expand the ho- 
mogeneous solutions to Eq. (4.301 to 0{q^). Define for 
convenience 



The solutions 0^*-' and -tpY' are given by 
,3 



(B7) 



6f) 



= 0, 



Ba) 
(B8b) 



32 



2 tanh" 



In 



32 



3 

J- 



! tanh 



2 In-" 2 - In — In /_ + In /+ In 
16 



21n/ln^ -4Li2^ 



3. Scalar Mode 



For the scalar mode it is sufficient to expand the ho- 



mogeneous solutions to Eq. (4.341 to 0{q). Doing so one 
obtains 

^(0) _ 4(2~-3r^) ^4 _ (l-3r^)K-f<(l-3r^)] ^ ^ 



9(l-r2)2 



9(l-r2)2 



= 0. 



and 



^0 



(0) _u^- (l-Sr'K 



(2-3r2)4 ' 



(B13a) 
(B13b) 



(B14a) 



(1) 



4w3(2-3r2) 



(4+ In 4)u^ - 4 {r^ In 64-^4- In 4) 
- [u^+ (l-3r2) 4] ln/| • (B14b) 
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The corresponding small momentum limit of the Wron- 
skian is 



,3L,4 



(l-r2)2(2-3r2)/ 
2ir (r2 In 8- In 2+2 



4(l-3r3) 



9< 



(B15) 



Substituting t he ab ove expansions into the definition of 
A^Q^'' in Eq. (4.701 and performing the radial integral, 
one finds 



Az; 



(4) 



3r2 + 



° 2^X3^13^ (l-3r2)u;, 
irq 9r2(l— r2)2 



(27r)(5(w - • g) 



Uh 



4(l-3r2 



(B16) 



Inserting AZ^'^'^ into Eq. (4.54a), we obtain 

ATo ^ ^ ^^^i^ {27t)S{oj -vq) 

-zr(l+z)2) r2(2-3r2+u2) 
{l-3r^)uhq ^ (l-3r2)2 
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